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PREFACE 

This  report  describes  the  work  performed  on  Air  Force  contract  F-40600 
-72-C-0015  by  the  Convair  Aerospace  Division  of  General  Dynamics  Cor¬ 
poration,  San  Diego  operation,  San  Diego,  California.  This  report  is 
identified  by  contractor's  number  CASD-AFS-73-008. 

The  work  was  administered  by  the  Department  of  the  Air  Force,  Head¬ 
quarters,  Arnold  Engineering  Development  Center  (TMP),  Arnold  Air 
Force  Station,  Tennessee.  Mr.  Ross  G.  Roepke,  AEDC  (DYX),  is  the 
Air  Force  technical  representative. 

This  program  was  conducted  in  the  research  and  engineering  department 
of  Convair  and  was  managed  by  S.  A.  Griffin  with  W.  K.  Alexander  as  co¬ 
ordination  engineer.  The  work  was  accomplished  between  April  and 
December  1973. 

Principal  contributors  to  the  study  include: 

Dynamics  R.  L.  Holt 

Aerodynamics  G.  J.  Fatta,  S.  P.  Tyler 

The  reproducibles  used  in  the  reproduction  of  this  report  were  supplied 
by  the  author. 
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SECTION  1 
INTRODUCTION 

The  international  concern  over  the  inability  of  existing  wind  tunnel  facilities  to  approach 
or  match  full-scale  Reynolds  numbers  has  led  to  the  development  of  a  variety  of  plans 
to  construct  high  Reynolds  number  transonic  wind  tunnels.  One  proposed  facility  is 
HIRT,  a  large  Ludwieg  tube  tunnel  with  an  8  by  10  foot  test  section.  The  facility  is  to 
be  located  at  Arnold  Engineering  Development  Center,  Arnold  Air  Force  Station, 
Tennessee  (References  1  and  2).*  The  High  Reynolds  Number  Transonic  Tunnel  (HIRT) 
will  operate  from  Reynolds  number  levels  of  present  wind  tunnels  to  full-scale  values 
for  future  aircraft. 

This  study  is  one  of  four  conducted  by  General  Dynamics  for  Arnold  Engineering 
Development  Center,  addressing  several  problem  areas  unique  to  the  proposed  HIRT 
facility.  An  initial  study  conducted  by  General  Dynamics  (Reference  3)  covered  the 
feasibility  of  designing  and  building  models  capable  of  withstanding  the  loads  and  en¬ 
vironmental  conditions  of  the  facility.  The  initial  study  concluded  that  models  capable 
of  being  tested  in  the  HIRT  facility  could  be  built  at  a  reasonable  cost  with  present-day 
techniques  and  materials.  The  additional  studies  have  examined  the  feasibility  of  multi¬ 
piece  flow-through  wind  tunnel  models  for  the  HIRT  facility,  (Reference  4),  six-com¬ 
ponent  internal  strain  gage  balance  designs  for  the  HIRT  facility  (Reference  5),  and  the 
expected  precision  of  test  data  in  the  HIRT  facility  (Reference  6). 


1.  J.  D.  Whitfield,  C.  J.  Schueler,  and  R.  F.  Starr,  ARO,  Inc. ,  "High  Reynolds 
Number  Transonic  Wind  Tunnels  —  Blowdown  or  Ludwieg  Tube",  AGARD-CP-83 
-71,  April  1971. 

2.  R.  G.  Roepke,  "The  High  Reynolds  Number  Transonic  Wind  Tunnel  HIRT  Proposal 
as  Part  of  the  National  Aeronautical  Facilities  Program,  "AIAA  Paper  No. 

72-1035,  September  1972. 

3.  W.  K.  Alexander,  S.  A.  Griffin  and  R.  L.  Holt,  General  Dynamics  Convair  Aero¬ 
space  Division,  "Wind  Tunnel  Model  Parametric  Study  for  Use  in  the  Proposed 

8  Ft  x  10  Ft  High  Reynolds  Number  Transonic  Wind  Tunnel  (HIRT)  at  Arnold 
Engineering  Development  Center,"  AEDC  Report  TR-73-47,  March  1973. 

4.  W.  K.  Alexander,  A.  E.  Brady  and  S.  A.  Griffin,  "Study  of  Multipiece  Flow-through 
Wind  Tunnel  Model  for  HIRT"  AEDC  Report  TR-75-60,  July  1974. 

5.  M.  L.  Kuszewski,  P.  J.  Mole  andS.  A.  Griffin,  "Study  of  Six-Component  Internal 
Strain  Gage  Balances  for  Use  in  the  HIRT  Facility,  "  AEDC  Report  TR-75-63, 

July  1974. 

6.  J.  Picklesimer,  and  W.  H.  Lowe,  "A  Study  of  Expected  Precision  In  the  Proposed 
AEDC  HIRT  Facility,"  AEDC  Report  TR-75-61,  July  1974. 

*  Since  completion  of  this  report  by  Convair,  a  final  decision  was  made  not  to  construct 
the  HIRT  at  AEDC  in  favor  of  a  continuous  cryogenic  wind  tunnel,  site  as  yet  un¬ 
determined. 
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The  high  load's  imposed  on  the  models  in  the  HIRT  facility  will  produce  considerable 
aeroelastic  deformations  to  the  primary  aerodynamic  surfaces.  An  undesired  wing 
twist  distribution  will  cause  a  disturbed  span  loading,  which  alters  the  configuration- 
induced  drag  characteristics.  Lift  and  pitching  moment  characteristics  are  also 
altered  and  lead  to  errors  in  trim  characteristics  caused  by  a  reduced  lift  curve  slope 
and  a  more  unstable  configuration. 

The  objective  of  this  task  was  to  study  the  effects  of  model  wing  deformations  obtained 
in  the  HIRT  facility.  These  deformations  were  analyzed  to  obtain  the  effect  on  the 
longitudinal  aerodynamic  characteristics  of  related  full  size  aircraft.  Two  study  con¬ 
figurations  were  analyzed.  First,  a  large  transport -type  aircraft,  the  Advanced 
Technology  Transport,  was  analyzed;  second,  a  variable-sweep  wing  supersonic  fighter 
type  aircraft  (F-lll)  was  analyzed. 

The  arrangement  of  this  report  is  as  follows.  Section  n  describes  the  study  program 
organization.  Section  HI  describes  the  test  conditions  that  were  considered  for  evalu¬ 
ating  the  aeroelastic  and  aerodynamic  characteristics.  Aeroelastic  methods  and 
analysis  are  presented  in  Section  IV.  For  the  aeroelastic  analysis,  a  modified  form 
of  the  Weissinger  L-method  was  used.  An  aerodynamic  influence  coefficient  matrix 
and  a  structural  flexibility  influence  coefficient  matrix  were  used  to  determine  the 
elastic  properties,  deformations,  and  rotations  of  the  airplane  structure.  Section  V 
provides  the  aerodynamic  methods  and  analysis,  with  a  description  of  the  study  air¬ 
craft.  Incremental  longitudinal  aerodynamic  characteristics  were  determined  from 
the  differential  aeroelastic  deformation  of  the  aircraft  and  the  wind  tunnel  model. 
Conclusions  are  presented  in  Section  VI  and  references  in  Section  VII.  An  additional 
study  was  conducted  to  obtain  the  aeroelastic  characteristics  of  an  airplane  composite 
wing  and  to  assess  the  HIRT  model/airplane  wing  deformation  differences  on  the 
longitudinal  aerodynamic  characteristics.  This  additional  study  is  presented  in 
Appendix  A. 
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SECTION  n 

PROGRAM  ORGANIZATION 

Convair  Aerospace  division  established  a  program  team  within  its  research  and  engi¬ 
neering  department  to  conduct  further  study  of  four  particular  problem  areas  associ¬ 
ated  with  testing  in  the  proposed  high  Reynolds  number  transonic  wind  tunnel  (HIRT). 
The  team  organized  to  conduct  the  studies  was  led  by  Mr.  S.  A.  Griffin,  who  reports 
to  Mr.  W.  T.  MacCarthy,  manager  of  wind  tunnels.  The  entire  organization  is  under 
the  senior  management  of  the  division  vice  president  of  research  and  engineering, 

Mr.  R.  H.  Widmer. 

The  program  operations  chart,  Figure  1,  illustrates  the  flow  of  information  from  the 
various  technical  groups  to  the  final  report. 


Figure  1.  Program  Organization 
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SECTION  m 
TEST  PLAN 


3. 1  TUNNEL  USE 

It  is  anticipated  that  a  great  demand  for  test  time  in  the  HERT  facility  will  result  in  a 
situation  where  testing  will  necessarily  be  restricted  to  tests  that  are  very  Reynolds 
number  dependent.  Hence,  use  of  the  HERT  facility  will  be  somewhat  unlike  that  of  cur¬ 
rent  transonic  wind  tunnels  in  that  plans  will  probably  be  much  more  operating  envelope 
oriented  or  will  concentrate  on  one  particular  test  problem.  It  is  current  practice  to 
test  a  configuration  without  particular  regard  to  its  envelope  so  that  data  is  obtained 
beyond  the  limits  that  the  flight  vehicle  can  achieve.  Such  practice  cannot  be  justified 
for  HERT  testing  because  of  the  impact  of  high  model  loads  and  model  deformations. 
HIRT  testing  will  best  be  carried  out  by  testing  in  critical  areas  of  the  operating 
envelope  to  obtain  the  configuration  performance  (typically  design  point  cruise  and 
maneuvering  parameters). 


3.  2  OPERATING  ENVELOPES 


For  the  purpose  of  this  study,  selected  points  within  the  operating  envelope  were  con¬ 
sidered  for  analysis  with  a  typical  cruise  condition  selected  to  establish  the  airplane 
wing  deflection  and  the  jig  twist  for  the  HIRT  models.  The  model  test  conditions  were 
considered  at  full-scale  Reynolds  number  with  the  tunnel  operating  nominally  at  ambient 
conditions.  Reduced  Reynolds  numbers  and  temperature  effects  were  also  considered. 
Reynolds  number  ranges  for  the  aircraft  covered  in  this  test  plan  are  illustrated  in 
Figure  2. 


Flight  vehicle  operational  envelopes  for  the  two  study  configurations  are  shown  in 
Figures  3  and  4.  Two  plots  are  presented  for  each  configuration.  The  first  plot  illus¬ 
trates  Rejy^Q^L  53  a  of  Mach  number  and  altitude  up  to  the  structural  limit 

load  factor.  The  second  plot  illustrates  the  HIRT  dynamic  pressures  for  the  equivalent 
envelope. 


3.  3  TEST  PLAN  DEVELOPMENT 

Tables  1  and  2  respectively  present  sum¬ 
maries  of  the  model  and  full-size  condi¬ 
tions  that  were  analyzed  to  determine : 

a.  Wing  distortions  (vertical  deflection 
and  twist). 

b„  Wing  bending  moments. 

c.  Wing  pitching  moments. 

d.  Wing  shear  loads. 

e.  Wing  loading. 

f.  Wing  section  properties  (El,  GJ). 

g.  Model  angle  of  attack. 


Figure  2.  Reynolds  Number  Ranges  for 
Aircraft  in  the  Test  Plan 
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Figure  3.  ATT  Operating  Envelope 


TUNNEL  FREESTREAM  DYNAMIC  PRESSURE  AT  300°K,q  (1000  lb /ft2) 


TUNNEL  FREESTREAM  DYNAMIC 
PRESSURE  AT  240°  K  (1,000  lb /ft2) 


Table  1.  ATT  HIRT  Model  Test  Plan 


AEDC-TR-75-62 


E>0 

O 


Tunnel  Wing 

Computer  Mach  Altitude  Load  Flight  Flight  Model  Tunnel  Tunnel  Tunnel  Storage  struct 
run  number  10s  factor  dyn  pres  c  a  dyn  pres  temp  R©/ft  Pressure  cross  sect 

no.  M  flt)  (g)  (psf)  L  (deg)  (psf)  fK)  (106)  (pBia)  (%) 


A1  0. 

90  30  1.0  357  0.366  4. 90  4248  300  30 

8  150  i( 

A2 

1.5 

0.549  7.20 

A3 

2.0 

0.732  9,40 

B1 

1 

0  1.0  8 

24  0.158  2.40  7891 

57. 7  260 

B2 

2.0 

0.316  4.4o 

B3 

3.0 

0.474  6.40 

B4 

4.0 

0.631  8.35 

B5 

5.0 

0.789  10.30 

Cl  0. 

70 

1.0  498  0.261  3.80  6256 

44 

9  330 

C2 

2.0 

0.522  7.70 

C3 

3.0 

0.783  11.50 

C4 

3.4 

0.887  13.20 

, 

D1  0. 

90  2 

0  1.0  551  0.236  3.40  5900 

43 

3  205 

D2 

2.0 

0.473  6.40 

D3 

3.0 

0.709  9.40 

D4 

4.0 

0.945  12.40 

'  ’ 

El  0. 

70 

1.0  333  0.390  5.50  3700 

31 

5  170  i 

E2 

1.5 

0.586  8.10 

E3 

2.0 

0.781  11.60 

E4 

2.5 

0.976  14. 30 

FI  0. 

90  1 

0  1.0  824  0.158  2.40  5846 

43 

2  2 

5 

F2 

2.0 

0.316  4.40 

F3 

3.0 

0.474  6.40 

F4 

4.0 

0.631  8.35 

F5 

5.0 

0.789  10.30 

G1 

1.0 

0.158  2.40  4100 

28. 

9  140 

G2 

2.0 

0.316  4.40 

G3 

3.0 

0.474  6.40 

G4 

4.0 

0.631  8.35 

. 

G5 

5.0 

0.789  10.30 

HI 

1.0 

0.158  2.40  58 

£6  240  57. 7  200 

H2 

2.0 

0. 316  4. 40 

H3 

3.0 

0.474  6.40 

H4 

4.0 

0.631  8.35 

HS 

5.0 

0. 789  10.  30 

Table  2.  F-lll  HIRT  Model  Test  Plan 


AEDC-TR-75-62 


AEDC-TR-75-62 


SECTION  IV 

AEROELASTIC  ANALYSIS 


4.1  METHODS 

Considerable  aeroelastic  deformation  of  models  in  the  HIRT  facility  is  a  certainty. 

A  means  of  accounting  for  the  model  aeroelasticity  is  necessary  before  the  test 
data  can  be  applied  to  the  full  size  design  problem.  A  simultaneous  solution  of 
equations  that  contain  both  the  aerodynamic  influence  functions  and  the  structural 
influence  functions  are  necessary  for  this  task.  A  method  based  on  a  modification 
of  the  Weissinger  L-Method  was  used  to  accomplish  this  analysis.  The  method  and 
digital  program  (Reference  7)  was  written  for  computing  steady-state  spanwise  load 
distributions  on  an  elastic  airplane  wing  for  specified  airplane  weights  and  load 
factors.  The  Weissinger  L-Method,  originally  developed  for  subsonic  flow,  is  valid 
for  supersonic  flow  providing  that  the  flow  over  the  wing  is  subsonic  (this  is  a  case 
where  the  leading  edge  of  a  swept  wing  is  behind  the  shock  cone). 

The  modified  Weissinger  L-Method  used  in  this  analysis  was  from  the  work  of  Gray 
and  Schenk  (Reference  8).  Additional  information  was  used  from  the  works  of 
Schindel  (Reference  9),  Borland  (Reference  10)  and  Blackwell  (Reference  13). 

The  aeroelastic  program  was  used  to  determine  the  loads  for  each  flight  condition 
of  the  airplane  and  the  resulting  wing  twist  and  bending.  Model  deformation  and 
loads  were  also  determined  for  the  equivalent  Reynolds  number  test  condition  in 
the  HIRT  facility.  Prior  to  the  determination  of  the  airplane/model  loading  conditions, 
rigid  aerodynamic  characteristics  of  the  airplane  configurations  were  determined 
through  the  use  of  wind  tunnel  test  data  and  DATCOM  (Reference  12). 


7.  "Application  of  Lifting  Line  Theory  to  Aircraft  Aeroelastic  Loads  Analysis," 
General  Dynamics  Convair  Report  No.  GDC-ERR- AN-1128,  February  1968.. 

8.  W.  L.  Gray  and  K.  M.  Schenk,  Boeing  Aircraft  Co.  ,  "A  Method  for  Calculating  the 
Subsonic  Steady-State  Loading  on  an  Airplane  with  a  Wing  of  Arbitrary  Planform 
and  Stiffness,"  NACA  TN  3030,  December  1953. 

9.  L.  H.  Schindel,  "An  Evaluation  of  Procedures  for  Calculating  Aerodynamic  Loads," 
AFFDL-TR-65-18,  May  1965. 

10.  C.  J.  Borland,  "Methods  of  Calculating  Aerodynamic  Loads  on  Aircraft  Structures," 
AFFDL-TR -66-37,  August  1966. 

11.  J.  A.  Blackwell,  Jr.  ,  "A  Finite  Step  Method  for  Calculation  of  Theoretical  Load 
Distribution  for  Arbitrary  Lifting  -  Surface  Arrangements  at  Subsonic  Speeds," 
NASA  TND-5335,  1969. 

12.  "USAF  Stability  and  Control  DATCOM,"  AFFDL. 
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Each  flight  condition  presented  in  Tables  1  and  2  was  analyzed  through  the  use  of  the 
aeroelastic  program.  This  section  presents  the  spanwise  variation  of  the  elastic 
wing  twist  and  the  wing  vertical  deflection  for  the  study  configurations.  Estimated 
incremental  aerodynamic  characteristics  are  presented  in  Section  V  of  this  report 
for  the  resulting  wing  twist  differences  between  the  airplane  and  the  HIRT  model 
under  equivalent  Reynolds  number  testing.  Model  and  full-size  geometric  data  used 
in  the  aeroelastic  analysis  are  presented  in  Figure  5  with  an  illustration  of  the  wing 
geometry  strip  arrangement. 
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Figure  5  Aerolastic  Program  Geometry 
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Figure  5  Aerolastic  Program  Geometry  (eont'd) 
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Figure  5  Aerolastic  Program  Geometry  (cont'd) 
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4.2  AEROELASTIC  EQUATIONS 


The  three  most  basic  equations  in  this  aeroelastic  analysis  are  as  follows. 


Kutta-Joukowsky  Relationship 


a  =  p  v  r 

where 

&  =  span  loading 
p  =  fluid  density 
Biot-Savart  Law  of  Vorticity 


V 

r 


=  velocity 


ds  =  circulation 


This  law  was  applied  to  the  bound  vortex  and  the  trailing  vortices  at  each  wing  strip  to 
calculate  the  aerodynamic  influence  functions: 

r(cos  a  -  cos  8) 

WP  4FR 

where  a  and  8  are  the  angles  between  the 
direction  of  the  vortex  segment  and  lines 
joining  the  ends  of  the  segment  to  the 
control  point  (sketch).  The  control  point 
in  this  theory  is  located  at  the  .75  e 
location  of  each  strip.  Wp  is  the  velocity 
induced  normal  to  the  plane  of  the  wing. 


r 


P 

PLAN  VIEW  END  VIEW 


Castigliano  Theorem 


A  part  of  the  general  Castigliano  Theorem  was  used  that  related  the  twisting  of  a 
structure  to  the  applied  bending  and  torsional  moments  by  means  of  virtually  applied 
moments.  Structural  properties  are  induced  by  the  El  and  GJ  terms. 


as  = 


f 


m  M  ds 


El 


+ 


t  T  ds 
GJ 


where 

Q!g  =  Elastic  angle  of  attack  change  due  to  bending  and  torsional  moments 
along  the  elastic  axis. 

m  =  Beam  bending  moment  per  unit  pitching  moment  applied  at  the  station 
at  which  0  is  to  be  determined. 

t  =  Torsional  moment  around  the  elastic  axis  per  unit  pitching  moment. 
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ds  =  Incremental  distance  along  elastic  axis. 

El  =  Effective  beam  stiffness  around  the  axis  of  the  bending  moments  M  and  m. 

GJ  =  Effective  torsional  stiffness  around  the  axis  of  the  torsional  moments 
T  and  t. 

Angle-of -attack  definitions  and  positive  sign  convention  are  as  follows^ 


✓  EQUILIBRIUM  POSITION  OF 
SECTION  ZERO-LIFT  LINE 


ROOT  SECTION  ZERO-LIFT  LINE 


Angle  of  Attack  at  Each  Strip: 


{aF}  =  (ar)  +  [a]  +  {ag} 


(l) 


where: 


[otF}  = 


[I]  -  QS  [EA]  + 


[AV„ 


LVZ 


4moj 


:  |j]A^]  +  [ApOjx  {cgc} 


r  i  P  mMds  p  tTds  _  r  . 

W-  /  ~E T~+  f  “gj"  =  q  [s]  x 

''a 


El  /  GJ 

~o  o 


£“gj = 


AV„ 


LVZ  J 


X  (“r )  -  H/PT}  X  PT  +  aT  +  (a0Lroot)  ~  K  Lsect} 


+  ^ncfyg)  +  ao 


'root 


a 


®SL 


6a  T  +  Oi  +  a 

SL  h  gn 


[AVz] 

a°Lbody  +  Qlinc 
sect  b°dy 

[  a  ] 

W 

[  a 

LVz  j 

a°Lroot+0!inc 

fus 

y  +  1 

15flI 

[6SP 

SP 


(ar)  =  root  angle  of  attack 
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4.3  TUNNEL  OPERATING  CONDITIONS 


Tunnel  operating  conditions  were  obtained  from  a  set  of  estimated  flow  property  curves 
for  the  proposed  HIRT  facility  (Reference  13).  These  curves  were  supplied  by  the  gov¬ 
ernment  for  use  in  this  study.  Tunnel  dynamic  pressures  as  a  function  of  Reynolds  num¬ 
ber  per  foot  are  presented  in  Figures  6  and  7  for  storage  temperatures  of  300°K  and 
240° K.  Since  Reynolds  number  is  sensitive  to  freestream  temperature,  a  significant 
reduction  of  dynamic  pressure  can  be  obtained  by  using  a  cooled  freestream  temperature. 


REYNOLDS  NO. /FT 


Figure  6.  Tunnel  Reynolds  Number 
and  Dynamic  Pressure  Characteristics 
at  300°K 


REYNOLDS  NO.  /FT 


Figure  7.  Tunnel  Reynolds  Number  and 
Dynamic  Pressure  Characteristics 
at  240°K 


13.  "Curves  of  Flow  Properties  for  HIRT  Operation,"  USAF,  June  1972. 


28 


AEDC-TR-75-62 


4.4  F-lll  AIRPLANE  AEROELASTIC  ANALYSIS 

The  F-lll  airplane  wing  deformation  characteristics  were  obtained  from  data  pre¬ 
sented  inReferences  14,  15  and  16.  Structural  properties  (El  and  GJ),  elastic  axis 
location,  airplane  geometry,  and  inertia  distribution  were  obtained  from  Reference  16. 
The  wing  bending  and  torsional  stiffness  properties  were  modified  to  account  for  the 
leading  and  trailing  edge  portions  of  the  wing  section.  The  wing  bending  moment  of 
inertia  and  the  pblar  moment  of  inertia  characteristics  of  the  F-lll  airplane  are  pre¬ 
sented  in  Figure  8  for  a  wing  sweep  of  50  degrees. 

Figure  9  shows  the  airplane  wing  jig  twist  distribution  for  the  50  degree  wing  sweep. 
Using  the  data  presented  in  Reference  16  and  the  data  shown  in  Figures  8  and  9,  the 
aeroelastic  program  was  used  to  obtain  the  deformation  characteristics  of  the  F-lll 
airplane  at  selected  flight  conditions. 

Figure  10  shows  the  elastic  wing  twist  distribution  for  the  F-lll  airplane  at  an 
altitude  of  10,000  feet  and  Mach  0.90  for  several  flight  load  factors.  Figure  11 
illustrates  the  elastic  wing  twist  distribution  as  a  function  of  load  factor  for  a  flight 
condition  of  20,  000  feet  and  Mach  0.90.  Figure  12  illustrates  the  elastic  wing  twist 
distribution  at  an  altitude  of  30,000  feet  and  Mach  0.90  for  several  flight  load  factors. 
From  Figures  10,  11  and  12  at  a  load  factor  of  2.0,  the  effect  of  increasing  dynamic 
pressure,  as  shown  in  Figure  13,  results  in  a  more  negative  wing  twist  (wash-out). 
Figure  14  illustrates  the  non-dimensional  vertical  deflection  of  the  airplane  wing  at 
load  factors  of  1.0,  2.0,  3.0  and  4.0  for  a  flight  condition  at  20, 000  feet  and  Mach 
0.90.  A  flight  condition  at  20,  000  feet  and  Mach  0.70  was  analyzed  resulting  in  the 
wing  elastic  distributions  shown  in  Figure  15.  Figure  16  shows  the  wing  elastic  dis¬ 
tributions  as  a  function  of  the  load  factor  for  the  airplane  at  an  altitude  of  10, 000  feet 
and  Mach  0.70. 


14.  "F-lll  Wing  Aeroelastic  Analysis  for  the  Safe  Life  Steel  Wing  Carry-Through- 
Box,"  General  Dynamics  Convair  Aerospace  Division,  Report  No.  FZS-12-1030, 
January  1970. 

15.  "F-111C  Wing  Aeroelastic  Analysis  for  the  Safe  Life  Steel  Wing  Carry-Through- 
Box, "  General  Dynamics  Convair  Aerospace  Division  Report  No.  FZS-12-5003, 
April  1971. 

16.  "F-lll  A/B  Basic  Data  Report,  "  General  Dynamics  Convair  Aerospace  Division 
Report  No.  FZS-12-157,  January  1966. 
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Figure  8.  Wing  Bending  Moment  of  Inertia  and  Polar  Moment 
of  Inertia  for  the  F-lll  Wing. 
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Figure  13.  Effect  of  Dynamic  Pressure  on  the  F-lll  Wing  Twist 
at  Mach  0.90  for  a  Load  Factor  of  2. 0 
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Figure  14.  Effect  of  Load  Factor  on  the  F-lll  Wing  Vertical  Deflection 
Characteristics  at  Mach  0.  90  (Altitude  =  20,  000  ft. ) 
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Figure  15.  Effect  of  Load  Factor  on  the  F-lll  Wing  Twist  at 
Mach  0.70  (Altitude  =20,000  ft. ) 


Figure  16.  Effect  of  Load  Factor  on  the  F-lll  Wing  Twist 
at  Mach  0.  70  (Altitude  =  10, 000  ft.) 


Table  3  gives  the  shear,  bending  moment,  and  pitching  torque  of  the  wing  at  the 
centerline  for  the  conditions  of  Table  2. 


Condition 

A 

B 

C 

D 

E 

No. 

Shear  Load, 

Air  +  Inertia  (thousands  of  lb) 

1 

23.  57 

30.28 

24.  72 

26.34 

22.27 

2 

32.  07 

46.50 

40.  87 

43.01 

30.48 

3 

40.  58 

62.  73 

57.01 

59.69 

38.  69 

4 

- 

78.95 

63.46 

76.36 

46.  90 

5 

- 

95.17 

- 

- 

Bending  Moment  (millions 

of  in-lb) 

1 

5.  50 

6.  73 

5.71 

6.02 

5.23 

2 

7.  56 

10.58 

9.61 

10.  03 

7.  24 

3 

9.  63 

14.44 

13.  50 

14.04 

9.23 

4 

- 

18.30 

15.  06 

18.05 

11.  23 

5 

- 

22. 15 

- 

“ 

- 

Pitching  Torque  (thousands  of  in-lb) 

1 

0.90 

*  0.98 

0.  90 

0.93 

0.  87 

2 

1.  29 

1.  72 

1.  64 

1.70 

1.  25 

3 

1.  68 

2.46 

2.39 

2.46 

1.  63 

4 

- 

3. 19 

2.  69 

3.22 

2.  01 

5 

- 

3.93 

- 

- 

- 

Table  3.  F-lll  Airplane  Shear,  Bending  Moment  and  Pitching  Torque 
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4.5  ATT  AIRPLANE  AEROELASTIC  ANALYSIS 

The  wing  bending  moment  of  inertia  and  the  polar  moment  of  inertia  characteristics  of 
the  ATT  airplane  are  presented  in  Figure  17.  Figure  18  shows  the  airplane  wing  jig 
twist  distribution  for  the  design  cruise  point.  This  distribution  produces  the  most 
favorable  span  loading  condition. 

Figure  19  shows  the  elastic  wing  twist  distribution  for  the  ATT  airplane  at  an  altitude 
of  36,000  feet  and  Mach  0.98  for  several  load  factors.  Figure  20  illustrates  the  elas¬ 
tic  wing  twist  distribution  as  a  function  of  load  factor  for  a  flight  condition  of  25,  000 
feet  and  Mach  0.80.  Figure  21  illustrates  the  non-dimensional  vertical  deflection  of 
the  airplane  wing  at  load  factors  of  1.0,  1.5,  2.0,  and  2.5  for  a  flight  condition  at 
25, 000  feet  and  Mach  0. 80.  The  vertical  deflection  is  interrelated  with  the  elastic 
wing  twist,  but  is  not  critical  to  the  wing  design  since  the  drag  increment  is  more 
readily  identified  with  the  wing  lift. 


Figure  17.  Wing  Bending  Moment  of  Inertia  and  Polar  Moment 
of  Inertia  for  the  ATT  Wing 
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Figure  18.  Wing  Jig  Twist  Distribution  for  the  ATT  Airplane 


Figure  19.  Effect  of  Load  Factor  on  the  ATT  Wing  Twist  at  Mach  0.98 
< Altitude  =  36,000  ft) 


Figure  20.  Effect  of  Load  Factor  on  the  ATT  Wing  Twist  at  Mach  0.80 
(Altitude  =  25,000  ft) 
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Figure  21.  Effect  of  Load  Factor  on  the  Wing  Vertical  Deflection 
of  the  ATT  Airplane  at  Mach  0. 80  (Altitude  =25,  000  ft) 


The  effect  of  load  factor  on  the  wing  twist  of  the  ATT  airplane  is  presented  in 
Figure  22  for  a  flight  condition  of  Mach  0. 60  at  an  altitude  of  10,  000  feet.  Figure  23 
illustrates  the  wing  twist  of  the  ATT  airplane  at  Mach  0.95  and  at  an  altitude  of 
20,000  feet.  Figure  24  shows  the  wing  twist  as  a  function  of  load  factor  for  a  flight 
condition  of  Mach  0.50  at  an  altitude  of  20,  000  feet  for  the  ATT  airplane.  The 
effect  of  dynamic  pressure  at  a  load  factor  of  2. 5  for  the  ATT  airplane  is  presented 
in  Figure  25.  Components  of  the  elastic  wing  twist  are  also  shown.  The  twist 
due  to  aerodynamic  loading  decreases  the  washout  with  increasing  dynamic  pressure 
since  the  center  of  pressure  is  forward  of  the  elastic  axis.  Wing  section  pitching 
moment  and  wing  inertia  lead  to  an  increase  in  the  wing  twist  as  the  dynamic 
pressure  increases.  The  resulting  twist  distribution  at  each  load  factor  becomes 
nearly  independent  of  dynamic  pressure.  This  is  more  of  a  coincidence  of  the 
combination  of  wing  sweep,  elastic  axis  location,  airfoil  characteristics,  and 
structural  properties  than  the  expected  result  of  a  vehicle  to  be  selected  at  random 
for  test. 
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Figure  25.  Effect  of  Dynamic  Pressure  on  the  ATT  Airplane  Wing  Twist 
at  a  Load  Factor  of  2. 5  (Altitude  *  20,  000  ft) 


The  shear,  bending  moment  and  torque  for  the  ATT  airplane  is  presented  in  Table 
4.  The  shear  loads  and  moments  are  determined  along  the  wing  elastic  axis. 


Table  4.  ATT  Airplane  Shear,  Bending  Moment  and  Pitching  Torque 


Condition 

A 

B 

C 

E 

G 

No.* 

Shear  Load, 

Air  +  Inertia  (thousands  of  lb) 

1 

112.  8 

127.5 

126.  0 

145.7 

112.  6 

2 

156. 1 

177.0 

174.4 

194.9 

161.9 

3 

199.5 

226.6 

222.  8 

244. 1 

211. 1 

4 

282.0 

276.1 

271.  2- 

293.3 

260.4 

Bending  Moment  (millions  of  lb) 

1 

60.46 

68.01 

68.27 

70.07 

64.  29 

2 

85.55 

95.90 

95.  74 

96.  52 

92.97 

3 

110.  64 

123.78 

123. 22 

122.  95 

121.  65 

4 

152.90 

151. 66 

150.  70 

149.4 

150.33 

Pitching 

Torque  (millions  of  in-lb) 

1 

-27.  73 

-29. 86 

-30.45 

-37. 92 

-23.  80 

2 

-35.  86 

-38.34 

-38.  89 

-46.  09 

-32.48 

3 

-43. 98 

-46. 83 

-47.  33 

-54. 25 

-41. 16 

4 

-54.  26 

-55.31 

-55. 77 

-62.41 

-49. 84 

*See  Table  1,  p.  23. 
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4.6  F-lll  HIRT  MODEL  AEROELASTIC  ANALYSIS 

The  wing  bending  moment  of  inertia  and  the  polar  moment  of  inertia  characteristics  of 
the  1/12  scale  F-lll  HIRT  model  are  illustrated  in  Figure  26.  Figure  27  illustrates 
the  model  jig  twist  distribution  which  causes  the  combination  of  elastic  twist  and  jig 
twist  of  the  model  to  equal  that  of  the  F-lll  airplane  at  the  design  cruise  point. 

The  effect  of  load  factor  on  the  F-lll  HIRT  model  wing  elastic  twist  and  vertical  de¬ 
flection  is  presented  in  Figures  28  through  35.  Effect  of  dynamic  pressure  on  the 
elastic  wing  twist  is  presented  in  Figure  36  for  the  F-lll  HIRT  model  at  Mach  0.90. 
The  effect  of  section  pitching  moment  (CmQ  )  is  less  than  the  aerodynamic  loading  due 
to  angle  of  attack  which  results  in  a  more  positive  wing  twist  with  increasing  dynamic 
pressure. 


Figure  26.  Wing  Bending  Moment  of  Inertia  and  Polar 
Moment  of  Inertia  for  the  F-lll  HIRT 
Model  Wing 


39 


ELASTIC  WING  TWIST  (deg)  ELASTIC  WING  TWIST  (deg)  MODEL  JIG  TWIST  (deg) 


AEDC-TR -75-62 


SEMISPAN 

Figure  28.  Effect  of  Load  Factor  on  the  F-lll  HIRT  Model  Wing  Twist 
at  Mach  0.90  and  Full  Scale  Reynolds  Number 
(Altitude  =  30,000  ft) 


Figure  29.  Effect  of  Load  Factor  on  the  F-lll  HIRT  Model  Wing  Twist 
at  Mach  0.90  and  Full  Scale  Reynolds  Number 
(Altitude  =  10,000  ft) 
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Figure  30.  Effect  of  Load  Factor  on  the  F-lll  HLRT  Model  at  Mach  0.70 
and  Full  Scale  Reynolds  Number  (Altitude  *  10,  000  ft) 
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Figure  31,  Effect  of  Load  Factor  on  the  F-lll  HIRT  Model  at  Mach  0.90 
and  Full  Scale  Reynolds  Number  (Altitude  =*  20, 000  ft) 


SEMISPAN 

Figure  32.  Effect  of  Load  Factor  on  the  Wing  Vertical  Deflection  of  the 
F-lll  HIRT  Model  at  Mach  0.90  and  Full  Scale  Reynolds 
Number  (Altitude  =  20,000  ft) 
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Figure  33.  Effect  of  Load  Factor  on  the  F-lll  HIRT  Model  Wing  Twist 
at  Mach  0.70  and  Full  Scale  Reynolds  Number 
(Altitude  =  20,000  ft) 


Figure  34.  Effect  of  Load  Factor  on  the  F-lll  HIRT  Model  Wing  Twist 


at  Mach  0.90  and  75%  Full  Scale  Reynolds  Number 
(Altitude  =  10,000  ft) 


SEMISPAN 

Figure  35.  Effect  of  Load  Factor  on  the  F-lll  HIRT  Model  Wing  Twist 
at  Mach  0.90  and  50%  Full  Scale  Reynolds  Number 
(Altitude  =  10,000  ft) 
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Figure  36.  Effect  of  Dynamic  Pressure  on  the  F-lll  HIRT  Model  Wing 
Twist  at  Mach  0.90  and  Full  Scale  Reynolds  Number 


A  relatively  large  effect  of  Reynolds  number  reduction  is  shown  in  Figure  37.  A  50 
percent  decrease  in  Reynolds  number  decreases  the  twist  by  40  percent. 

Table  5  presents  the  shear,  bending  moment,  and  pitching  torque  of  the  F-lll  1/12 
scale  model  at  the  wing  centerline  for  the  conditions  of  Table  2. 


SEMISPAN 


Figure  37.  Effect  of  Reynolds  Number  on  the  F-lll  HIRT  Model  Wing 
Twist  at  Mach  0.90  and  at  a  Load  Factor  of  5.0 
(Altitude  =  10,000  ft) 
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Table  5.  F-lll  HIRT  Model  Shear,  Bending  Moment  and  Pitching  Torque 


Condition 

A 

B 

C 

D 

E 

F 

G 

No. 

Shear  Load,  Air  +  Inertia  (thousands  of  lb) 

1 

2.  24 

2.32 

2. 13 

2.28 

2.06 

L  78 

1.28 

2 

3.  08 

3.64 

3.62 

3.  80 

2.86 

2.  80 

2.03 

3 

3.  93 

4.98 

5. 10 

5.33 

3.65 

3.82 

2.78 

4 

- 

6.30 

5.  69 

6.  86 

4.45 

4.  85 

3.  52 

5 

- 

7.  63 

- 

- 

- 

5.  88 

4.27 

Bending  Moment  (thousands  of  in-lb) 

1 

43.  92 

43.74 

41.58 

43.  89 

40.67 

33.  82 

24.  59 

2 

61. 12 

70. 35 

71.60 

74.  69 

56.96 

54.50 

39.74 

3 

78.38 

97.07 

101.  63 

105.56 

73.24 

75.26 

54.  94 

4 

- 

123.  69 

113.  61 

136. 44 

89.53 

95.  94 

70.  09 

5 

- 

150.31 

- 

- 

- 

116.  62 

85.  24 

Pitching  Torque  (thousands  of  in-lb) 

1 

6.  35 

5.78 

5.  83 

6. 10 

5.92 

4.  52 

3.  32 

2 

9.  08 

10.00 

10.58 

10.98 

8.50 

7.79 

5.  71 

3 

11.  81 

14.24 

15.34 

15.  87 

11.07 

11.  08 

8.  12 

4 

- 

18.46 

17.  23 

20.76 

13.65 

14.35 

10.  52 

5 

- 

22.68 

“ 

- 

17.  63 

12.91 

A  comparison  of  the  stiffness  ratio  (El/GJ)  for  the  F-lll  airplane  and  the  1/12  scale 
HIRT  model  is  presented  in  Figure  38. 


It  is  anticipated  that  the  stiffness  ratio  for  wings  fabricated  of  composite  materials 
will  approach  a  value  of  3.0.  The  aeroelastic  analysis  for  this  condition  is  pre¬ 
sented  in  Appendix  A  for  the  F-lll  wing. 


Figure  38.  Stiffness  Ratio  Characteristics  of  the  F-lll  Airplane  and 
HIRT  Model 
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The  effect  of  tunnel  storage  temperature  conditions  on  the  elastic  wing  twist  of  the 
F-lll  HIRT  model  is  presented  in  Figure  39  for  equivalent  Reynolds  number  test 
conditions.  The  reduction  in  this  temperature  from  300°  K  to  240°  K  results  in  about 
a  21  percent  reduction  in  the  wing  tip  twist  at  a  constant  load  factor.  This,  a  reduction 
in  tunnel  temperature  is  an  attractive  method  of  controlling  the  tunnel  dynamic 
pressure  and  therefore  the  model  wing  deformations  while  maintaining  equivalent  full 
scale  Reynolds  number. 
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Figure  39.  Effect  of  Tunnel  Temperature  on  the  F-lll  HIRT  Model  Wing 
Twist  at  Full  Scale  Reynolds  Number 


A  comparison  of  the  F-lll  airplane  and  the  HIRT  model  wing  vertical  deflection  is 
presented  in  Figure  40  for  an  equivalent  Reynolds  number  condition.  The  nondimen- 
sional  vertical  deflection  ratio  at  the  wing  tip  for  the  HIRT  model  is  about  18  percent 
greater  than  the  F-lll  airplane  at  a  constant  load  factor  and  an  equivalent  Reynolds 
number.  The  wing  vertical  deflection,  measured  along  the  elastic  axis,  has  been 
nondimensionalized  with  respect  to  the  wing  mean  aerodynamic  chord.  Wing  vertical 
deflection  is  an  integral  portion  of  the  resulting  wing  twist  distribution.  For  this 
case,  the  model  wing  deflection  and  twist  are  greater  than  the  airplane  at  an  equivalent 
Reynolds  number. 
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Figure  40.  Comparison  of  the  F-lll  Airplane  and  the  F-lll  HIRT 
Model  Wing  Vertical  Deflection  at  Equivalent  Reynolds 
Number. 
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4.  7  ATT  HTR.T  MODEL  AEROELASTIC  ANALYSIS 

The  wing  bending  moment  of  inertia  and  the  polar  moment  of  inertia  characteristics 
of  the  steel  ATT  HIRT  model  are  shown  in  Figure  41.  Figure  42  shows  the  model 
wing  jig  twist  distribution  which  causes  the  combination  of  the  elastic  twist  and 
jig  twist  of  the  model  to  equal  that  of  the  airplane  wing  twist  at  the  design  cruise 
point. 

Figure  43  shows  the  elastic  wing  twist  distribution  for  the  ATT  HIRT  model  at 
Mach  0.98  for  several  flight  load  factors.  Tunnel  conditions  of  dynamic  pressure 
and  temperature  produced  the  equivalent  full  scale  Reynolds  number  for  a  flight 
condition  of  Mach  0.98  at  an  altitude  of  36,  000  feet. 

Figure  44  illustrates  the  elastic  wing  twist  distribution  for  the  ATT  HIRT  model 
as  a  function  of  load  factor  for  a  flight  condition  of  Mach  0. 80. 

Effect  of  load  factor  on  the  ATT  HIRT  model  for  a  flight  condition  of  Mach  0.60  is 
presented  in  Figure  45  for  full  scale  Reynolds  number  conditions.  Figures  46,47 
and  48  illustrates  the  effect  of  load  factor  on  the  ATT  HIRT  model  for  equivalent 
Reynolds  number  conditions . 

Effect  of  wing  stiffness  ratio  for  the  ATT  HIRT  model  wing  is  presented  in  Figure 
49.  A  solid  steel  wing  would  represent  a  100  percent  structural  arrangement. 

The  model  wing  has  been  designed  with  pressure  instrumentation  within  the  wing. 
This  results  in  a  structural  degradation  of  the  wing  cross-sectional  properties. 

The  basic  model  arrangement  has  a  65  percent  wing  structural  effectiveness.  For 
a  reduction  in  the  wing  structural  characteristics  to  52  percent  of  the  solid  cross- 
section  causes  the  wing  twist  to  increase  by  15  percent  over  the  basic  configuration 
as  shown  in  Figure  49. 

The  effect  of  a  tunnel  dynamic  pressure  at  a  load  factor  of  2.  5  for  the  ATT  HIRT 
model  is  presented  in  Figure  50.  Components  of  the  elastic  wing  twist  are  also 
shown.  The  twist  due  to  aerodynamic  loading  decreases  the  washout  with  increasing 
dynamic  pressure  since  the  center  of  pressure  is  forward  of  the  elastic  axis.  Wing 
section  pitching  moment  and  wing  inertia  lead  to  an  increase  in  the  washout. 

The  effect  of  tunnel  dynamic  pressure  on  the  ATT  HIRT  model  wing  vertical  deflec¬ 
tion  is  presented  in  Figure  51  for  a  load  factor  of  2.  5  at  the  equivalent  full  scale 
Reynolds  number.  For  the  constant  load  factor,  the  lower  dynamic  pressure 
condition  results  in  larger  vertical  deflection.  Wing  vertical  deflection  as  a  function 
of  the  load  factor  is  presented  in  Figure  52  for  the  ATT  HIRT  model  at  an  equivalent 
full  scale  Reynolds  number  condition. 
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Figure  43.  Effect  of  Load  Factor  on  the  ATT  HIRT  Model  at  Mach  0.98 
and  Full  Scale  Reynolds  Number  (Altitude  -36,000  ft) 


Figure  44.  Effect  of  Load  Factor  on  the  ATT  HIRT  Model  at  Mach  0.80  and  Full 
Scale  Reynolds  Number  (Altitude  =25,000  ft) 


Figure  45.  Effect  of  Load  Factor  on  the  ATT  HIRT  Model  at  Mach  0. 60 
and  Full  Scale  Reynolds  Number  (Altitude=10,  000  ft) 
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Figure  46. 


Effect  of  Load  Factor  on  the  ATT  HIRT  Model  at  Mach  0.98 
and  Full  Scale  Reynolds  Number  (Altitude® 36, 000  ft) 


SEMISPAN 


Figure  47.  Effect  of  Load  Factor  on  the  ATT  HIRT  Model  at  Mach  0.95 
and  Full  Scale  Reynolds  Number  (Altitude® 20, 000  ft) 


SEMISPAN 

Figure  48.  Effect  of  Load  Factor  on  the  ATT  HIRT  Model  at  Mach  0.50 
and  Full  Scale  Reynolds  Number  (Altitude® 20, 000  ft) 
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Figure  49. 


Effect  of  Wing  Stiffness  Ratio  for  the  ATT  HIRT  Model  Wing. 


Figure  50.  Effect  of  Dynamic  Pressure  on  the  ATT  HIRT  Model 
Wing  Twist  at  a  Load  Factor  of  2. 5. 
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Figure  51.  Effect  of  Tunnel  Dynamic  Pressure  on  the  ATT  HIRT  Model  Wing 
Vertical  Deflection  at  Full  Scale  Reynolds  Number. 


SEMISPAN 

Figure  52.  Effect  of  Load  Factor  on  the  ATT  HIRT  Model  Wing  Vertical  Deflection 
at  Mach  0.80  and  Full  Scale  Reynolds  Number  (Altitude®25, 000  ft). 
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Figure  53  compares  the  structural  stiffness  ratios  of  the  ATT  airplane  and  the  ATT 
HIRT  model.  The  ATT  HIRT  model  wing  stiffness  ratio  is  about  one  half  of  the 
ATT  airplane  wing. 


Figure  53.  Stiffness  Ratio  Characteristics  of  the  ATT  Airplane  and 
the  ATT  HIRT  Model. 


The  elastic  wing  twist  for  the  ATT  HIRT  model  is  compared  with  the  ATT  airplane 
elastic  wing  twist  in  Figure  54  for  the  design  cruise  point.  The  comparison  at 
equivalent  Reynolds  number  conditions  indicates  that  the  model  elastic  twist 
distribution  may  be  held  within-  0. 1  degrees  angle  of  attack  at  each  wing  section 
relative  to  the  airplane  wing. 
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Figure  54.  Comparison  of  the  ATT  Airplane  and  the  ATT  HIRT 
Model  at  the  Equivalent  Reynolds  Number  for  the 
Design  Cruise  Point 
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Table  6  shows  the  shear  load,  beading  moment  and  pitching  torque  for  the  ATT  HlftT 
model  along  the  wing  elastic  axis.  The  flight  conditions  are  presented  in  Table  1. 


Table  6.  ATT  HIRT  Model  Shear,  Bending  Moment  and  Pitching  Torque, 


Condition 

A 

B 

C 

D 

E 

G 

No. 

Shear  Load,  Air  +  Inertia  (thousands  of  lb) 

1 

5. 44 

5.  89 

5.48 

5.32 

6.  06 

5.30 

2 

7.  62 

8.  24 

7.66 

7. 45 

8.  15 

7.67 

3 

9.81 

10.  62 

9.  85 

9.  59 

10,26 

10.  04 

4 

13.65 

12.99 

12.  02 

13.  34 

12.35 

12. 42 

Bending  Moment  (thousands  of  in-lb) 

1 

96.35 

106.32 

100. 86 

92.48 

103.58 

100. 02 

2 

135. 00 

148. 43 

140.  09 

129. 62 

139.  29 

144. 23 

3 

173.  87 

190.  74 

179. 50 

166.  96 

175.31 

188.  52 

4 

241.  88 

233. 06 

218.55 

232.30 

211.  02 

232.  72 

Pitching  Torque  (thousands  of  in- 

-lb) 

1 

-29. 99 

-32.  61 

-31.61 

-29. 33 

-41.  42 

-22. 43 

2 

-34. 72 

-37.  79 

-36.48 

-33. 81 

-45. 64 

-28.07 

3 

-39.48 

-43. 01 

-41.37 

-38.31 

-49. 90 

-33. 73 

4 

-47.  81 

-48.  22 

-46.  22 

-46.  18 

-54. 12 

-39.38 

The  effect  of  tunnel  storage  temperature  conditions  on  the  elastic  wing  twist  of  the 
ATT  HIRT  model  is  presented  in  Figure  55  for  equivalent  Reynolds  number  test  con¬ 
ditions.  The  reduction  in  tunnel  temperature  from  300°  K  to  240°  K  results  in  about  a 
20  percent  reduction  in  the  wing  tip  twist  at  a  constant  load  factor.  Again  the  effect 
of  storage  temperature  as  a  possible  control  of  model  deformation  should  be  noted. 
Reduction  in  the  tunnel  temperature  may  best  be  applied  to  the  situation  where  the 
steel  model  is  structurally  limited  in  comparison  to  the  airplane  which  has  been  de¬ 
signed  for  high  load  factors. 
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Figure  55.  Effect  of  Tunnel  Temperature  on  the  ATT  HIRT  Model  Wing 
Twist  at  Full  Scale  Reynolds  Number. 


A  comparison  of  the  ATT  airplane  and  the  HIRT  model  wing  vertical  deflection  is  pre¬ 
sented  in  Figure  56  for  an  equivalent  Reynolds  number  condition.  The  nondimensional 
vertical  deflection  ratio  at  the  wing  tip  for  the  HIRT  model  is  about  31  percent  less 
than  the  ATT  airplane  at  a  constant  load  factor  and  an  equivalent  Reynolds  number. 
The  wing  deflection  is  an  integral  part  of  the  wing  twist  distribution.  For  this  case, 
the  model  wing  vertical  deflection  and  twist  are  less  than  the  airplane  at  an  equivalent 
Reynolds  number. 
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Figure  56.  Comparison  of  the  ATT  Airplane  and  the  ATT  HIRT  Model 
Wing  Vertical  Deflection  at  Equivalent  Reynolds  Number 
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4.  8  ITERATIVE  ANALYSIS  OF  WIND  TUNNEL  RESULTS 

The  correlation  of  the  full  size  analysis  with  the  scale  model  results  will  be  dependent 
upon  getting  the  model  to  deform  like  the  full  size  vehicle  as  well  as  maintaining  the 
Reynolds  number.  In  a  real-time  testing  situation,  it  may  be  necessary  to  change  the 
tunnel  dynamic  pressure  and  temperature  to  achieve  the  desired  deformation  for  a 
given  angle  of  attack  of  the  wing  root.  To  reconcile  the  differences  between  the  theo¬ 
retical  analysis  of  the  model  and  the  observed  HIRT  results,  it  will  be  desirable  to 
iteratively  analyze  the  effects  of  the  aerodynamic  spanwise  distributions  of  ,  Cm  , 
and 9  Cmcy4^  C^;  the  structural  spanwise  distributions  of  El  and  GJ;  and,  perhaps, 
the  spanwise  variation  of  the  elastic  axis.  As  the  test  series  progresses,  analyses 
can  identify  corrections  to  the  above  variables  required  to  match  calculations  with 
the  model  deformations. 

Matching  the  theoretical  model  analysis  to  the  deformations  measured  in  the  HIRT 
tunnel  is  important  because  the  aerodynamic  distributions  of  both  the  model  and  full 
size  vehicle  should  be  in  agreement.  The  structural  characteristics  of  a  full  size 
vehicle  are  normally  reliable  for  these  purposes  when  the  entire  structure  is  con- 
s  idered  as  opposed  to  using  the  wing  box  alone.  Should  the  full  size  wing  aerodynamics 
not  be  in  agreement  with  the  model,  the  deformations  will  be  erroneous.  Because  of 
the  need  to  analyze  and  test  the  model  and  full  size  vehicle  simultaneously,  it  is  desir¬ 
able  to  perform  the  analysis  rapidly. 

This  requirement  is  a  natural  application  for  the  interactive  graphics  cathode-ray 
tube  display.  The  use  of  a  real-time  mode  rather  than  a  conventional  batch  mode 
means  that  a  satisfactory  matching  of  the  aerodynamic  and  structural  parameters  and 
deformations  can  be  accomplished  in  minutes  while  the  test  article  is  in  the  tunnel, 
rather  than  by  days  of  analysis  interspersed  with  tunnel  setup  and  teardown  activity. 

The  cost  of  HIRT  tunnel  activity  will  be  much  greater  than  the  digital  analysis.  The 
latter  will  be  the  most  useful  when  operated  in  the  interactive  graphics  mode. 
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SECTION  V 

AERODYNAMIC  ANALYSIS 


The  high  loads  imposed  on  the  models  in  the  HIRT  facility  will  produce  aeroelastic 
deformation  that  will  be  different  from  those  of  the  flight  vehicle  due  to  differences 
in  the  material  characteristics  from  which  the  aircraft  and  models  are  constructed. 

The  procedure  used  during  this  study  considered  matching  the  model  wing  twist  with 
the  flight  vehicle  wing  twist  for  a  selected  point,  such  as  a  cruising  altitude  and  Mach 
number.  Additional  points  within  the  flight  envelope  were  analyzed  and  the  aeroelastic 
differences  determined. 

Rigid  aerodynamic  characteristics  were  determined  for  the  model  and  the  airplane 
for  each  deformed  wing  shape.  The  aerodynamic  methods  did  not  consider  the  distor¬ 
tion  components  due  to  tunnel  constraints  for  the  model.  The  linearized  aerodynamic 
methods  used  in  the  analysis  considered  the  wing  twist  as  a  linear  twist  distribution 
from  a  selected  point  on  the  wing  semispan  to  the  wing  tip.  Incremental  linearized 
aerodynamic  coefficients  (AC^.ACp  andACm)  were  then  determined  for  the  study 
configurations  at  each  of  the  selected  flight  conditions. 

5. 1  DESCRIPTION  OF  AIRCRAFT  AND  MODELS 

Two  types  of  aircraft  have  been  selected  to  be  analyzed  during  this  study.  The  aircraft 
represent  typical  type  vehicles  that  would  be  likely  candidates  for  testing  in  the  HIRT 
facility.  The  aircraft  configurations  were: 

a.  The  Advanced  Technology  Transport  (ATT),  which  is  a  highly  loaded,  large, 
swept-wing  transport  aircraft. 

b.  The  Convair  Aerospace  F-lll  fighter,  which  is  a  swing  wing  supersonic  combat 
aircraft. 

The  models  designed  for  the  HIRT  facility  were  based  on  a  solid-steel  construction  for 
sting  mounting. 

5. 1. 1  Advanced  Technology  Transport  (ATT) 

The  ATT  aircraft  used  in  this  study  was  a  predesign  configuration.  The  aircraft  was 
a  three-jet,  low-wing  transport  with  a  44-degree  swept-back  wing  using  a  super¬ 
critical  airfoil.  The  wing  had  an  aspect  ratio  of  8.  0,  a  theoretical  wing  area  of  2282 
square  feet  and  a  taper  ratio  of  0.390.  The  fuselage  was  area  ruled,  with  the  wing 
empennage  and  engines  arranged  to  achieve  the  optimum  cross -sectiontal  area  distri¬ 
bution  for  the  aircraft  cruising  at  a  design  Mach  number  of  0.  98.  Overall  dimensions 
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of  the  aircraft  included  the  length,  183.3  feet;  height,  51  feet;  and  span,  135. 1  feet. 
Dimensional  characteristics  of  the  ATT  are  presented  in  Table  7;  a  general  arrange¬ 
ment  is  presented  in  Figure  57. 

Table  7.  Dimensional  Characteristics  of  the  Advanced  Technology  Transport 


Wing 

Area  2282  ft2 

Span  135.1ft 

Root  Chord  24.162  ft 

Tip  Chord  9.416  ft 

Aspect  Ratio  8.  0 

Mean  Aerodynamic  Chord  20. 1  ft 

Leading  Edge  Sweep  44.35° 

Dihedral  3° 

Incidence  of  the  Root  0 

Taper  Ratio  0.3897 

Horizontal  Stabilizer 

.  ^ 

Area  395  ft 

Span  46.58  ft 

Root  Chord  16.33  ft 

Tip  Chord  6.58  ft 

Aspect  Ratio  3.00 

Mean  Aerodynamic  Chord  14. 1  ft 

Leading  Edge  Sweep  45° 

Dihedral  7° 

Taper  Ratio  0.4065 

Vertical  Stabilizer 

2 

Area  215  ft 

Span  18.  5  ft 

Root  Chord  17.  0  ft 

Tip  Chord  6.25  ft 

Aspect  Ratio  1.59 

Mean  Aerodynamic  Chord  24.36  ft 

Leading  Edge  Sweep  49.  5° 

Taper  Ratio  0,6992 


The  ATT  model  designed  for  the  HIRT  facility  (Reference  3)  was  a  l/24-scale  all 
sted  model.  The  general  arrangement  of  the  model  is  presented  in  Figure  58. 

The  wings  shown  in  Figures  57  and  58  are  similar  except  for  size.  The  model  wing 
had  allowances  for  pressure  instrumentation,  which  resulted  in  a  35  percent 
degradation  in  the  structural  strength  of  the  solid  steel  wing. 
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5.1.2  Variable  Swept  Wing  Fighter  (F-lll) 

The  F-lll  aircraft  used  in  this  study  was  a  production  design  F-111A  airplane.  The 
airplane  was  a  twin  jet  high  wing  combat  aircraft  with  a  variable  swept  wing  utilizing 
a  64A  series  airfoil.  The  aeroelastic  analysis  was  performed  with  the  wing  swept  to 
50  degrees.  The  wing  has  an  aspect  ratio  of  4.325,  a  taper  ratio  of  0.325  and  a  wing 
area  of  537.  5  square  feet  at  the  50-degree  leading  edge  sweep  position.  Overall  dimen¬ 
sions  of  the  aircraft  included  the  length,  75.54  feet;  height,  17.5  feet;  and  span,  48.22 
feet  at  the  50-degree  leading  edge  sweep.  Dimensional  characteristics  of  the  airplane 
are  presented  in  Table  8;  a  general  arrangement  in  Figure  9. 


Table  8.  Dimensional  Characteristics  of  the  Convair  Aerospace  F-lll  Aircraft 


Wing  at  50°  L.E. 

Sweep 

Area 

537.5  ft 

Span 

48.  22  ft 

Root  Chord 

12.56  ft 

Tip  Chord 

4.07  ft 

Aspect  Ratio 

4.325 

Mean  Aerodynamic  Chord 

9.04  ft 

Leading  Edge  Sweep 

O 

O 

LO 

Dihedral 

0“ 

Incidence  at  the  Root 

1— 1 

O 

o 

Taper  Ratio 

0.325 

Horizontal  Stabilizer 

Vertical  Stabilizer 

Area  154.  2  ft 

Area 

111.  7  ft2 

Span 

29.33  ft 

Span 

8.90  ft 

Root  Chord 

15.008  ft 

Root  Chord 

17.79  ft 

Tip  Chord 

5.012  ft 

Tip  Chord 

7.31  ft 

Aspect  Ratio 

2.07 

Aspect  Ratio 

0.709 

Mean  Aerodynamic  Chord 

11.  46  ft 

Mean  Aerodynamic  Chord 

13.  28  ft 

Leading  Edge  Sweep 

57.5° 

Leading  Edge  Sweep 

55° 

Dihedral 

0° 

Taper  Ratio 

0.411 

Taper  Ratio 

0.334 

The  F-lll  model  that  was  designed  for  the  HIRT  facility  was  a  l/12-scale  all-steel 
model.  The  wing  was  designed  as  a  nonpivoting  all-steel  structure  in  Reference  3.  A 
general  arrangement  of  the  F-lll  model  is  presented  in  Figure  60. 
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5. 2  AERODYNAMIC  METHODS 

Several  methods  have  been  developed  to  provide  a  quick  response  capability  for  evalu¬ 
ating  an  aircraft  design  and  performing  some  basic  trade  studies.  The  methods  used 
for  this  study  employed  existing  techniques  that  are  applicable  to  perform  a  data  corre¬ 
lation  study.  A  digital  computer  program  designated  the  Aeromodule  Procedure  (Refer¬ 
ence  17)  has  been  used  to  evaluate  the  lift  and  drag  characteristics  of  the  selected  study 
aircraft.  The  program  features  a  minimum  geometry  input,  an  airfoil  section  library, 
several  methods  for  minimum  drag  prediction  and  computes  a  lift  curve  slope,  an  angle 
of  attack,  minimum  drag,  drag  due  to  lift,  wave  drag  and  transonic  drag  rise. 

Wing  lift  curve  slope  was  predicted  by  use  of  an  empirical  approach  developed  to  in¬ 
clude  continuous  subsonic,  transonic  and  supersonic  application,  provide  realistic 
peak  lift-curve  slope  in  the  transonic  region,  approximate  transonic  bucket  trends, 
and  provide  first-order  wing-body  interference  effects.  The  total  lift  curve  slope  in¬ 
cluded  the  contributions  from  the  wing,  body,  and  tail.  The  wing  lift  curve  slope 
basically  started  with  the  Polhamus  equation  for  trapezoidal  wings: 
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The  incremental  term  ina^  due  to  twist  was  extrapolated  from  the  results  of  a  para¬ 
metric  study  reported  by  Gilman  and  Burdges  (Reference  18).  Their  results,  which 
neglect  body  effects,  are  valid  for  wings  with  linear-element  twist. 


17.  "Computer  Program  Details  of  Aeromodule  Version  III:  Lift  and  Drag  Analysis 
of  Aerodynamic  Configurations,"  General  Dynamics  Convair  Aerospace  Division 
Report  No.  MR-A-2089,  December  1970. 

18.  B.  G.  Gilman  and  K.  P.  Burdges,  "Rapid  Estimation  of  Wing  Aerodynamic 
Characteristics  for  Minimum  Induced  Drag,"  Journal  of  Aircraft,  Vol.  4, 

No.  6,  November  1969. 
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Drag  due  to  lift  was  determined  by  several  techniques,  depending  on  the  conditions  at 
which  the  solution  was  required.  For  subsonic  Mach  numbers,  drag  due  to  lift  was 
calculated  in  the  low-lift  parabolic  region,  in  the  nonlinear  stall  region,  and  in  the 
high  lift  region  above  stall.  The  drag  due  to  lift  term  is  basically  expressed  by 


k(cl-acl>2. 


The  polar  shape  factor,  K,  was  determined  from  the  wing  planform  geometry,  the 
span  efficiency  factor,  and  the  wing  leading  edge  suction  factor,  which  was  a  function 
of  the  leading  edge  radius,  Reynolds  number,  and  Mach  number.  The  polar  displace¬ 
ment  term,  ACL,  is  a  function  of  the  wing  design  lift  coefficient,  the  Mach  number, 
and  the  span  efficiency  factor. 

Using  the  output  from  the  Weissinger-L  method,  which  produced  lift,  pitching  moments 
and  aeroelastic  deformation  of  the  aircraft,  the  rigid  aerodynamic  program  was  used 
to  determine  the  drag  characteristics  of  the  aircraft  at  a  single  deformed  point.  By 
determining  the  rigid  aerodynamic  characteristics  for  various  wing  deformations,  the 
flexible  aerodynamic  characteristics  were  determined.  This  method  has  been  used  to 
evaluate  the  characteristics  of  the  flight  vehicle  and  the  model  with  the  differences 
being  presented  and  discussed  in  this  study. 


5. 3  LONGITUDINAL  AEROELASTIC/AERODYNAMIC  CHARACTERISTICS  OF  THE  ATT 

Aeroelastic  deformations  of  the  ATT  wing  were  analyzed  for  several  flight  conditions 
of  the  aircraft  and  the  corresponding  conditions  of  the  HLRT  model  through  the  use  of 
the  modified  Weissinger-L  method  (see  Chapter  4).  For  the  ATT  model,  the  deflected 
wing  of  the  aircraft  in  level  flight  (n=l)  at  36,000  feet  and  Mach  0.  98  was  assumed  to 
be  the  design  twist  used  to  determine  the  model  wing  jig  twist.  The  resulting  jig  twist 
as  shown  in  Figure  42  resulted  in  matching  the  model  twist  distribution  to  that  of  the 
aircraft  at  the  selected  condition  of  Mach  number  and  lift  coefficient.  The  model  wing 
with  the  thin  supercritical  airfoil  section  was  designed  to  be  constructed  of  steel  with 
35  per  cent  of  the  wing  section  not  considered  as  contributing  to  the  bending  and  tor¬ 
sional  stiffness  characteristics  because  of  high  lift  device  considerations  or  the  inclu¬ 
sion  of  pressure-tube  routing  channels. 

5.  3.  1  Characteristics  at  Mach  0.  98  and  36,000  Feet 

Figure  61  illustrates  the  comparison  of  the  full-size  wing  twist  with  that  of  the  model 
for  HIRT  test  conditions  simulating  full-scale  Reynolds  number  at  36,000  feet  and 
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Figure  61.  Comparison  of  the  ATT  Model  Wing 
Twist  with  the  ATT  Airplane  Wing 
Twist  for  Level  Flight  at  Mach  0.  98 


Mach  0.  98.  As  the  load  conditions  are  varied  by  pitching  from  the  wing  twist  matching 
point,  the  model  wing  deflections  diverge  from  the  airplane  wing  deflection.  Figure  62 
illustrates  the  comparison  of  the  airplane  wing  tip  twist  with  that  of  the  model  for  in¬ 
creasing  load  factors.  The  data  in  Figure  62  show  that  the  airplane  wing  twists  1.  7 
degrees  per  g,  while  the  model  wing  twists  1.  5  degrees  per  g  under  the  HERT  test  con¬ 
ditions,  The  data  indicates  that  the  model  is  staffer  under  the  IIIRT  test  conditions 
than  the  airplane  under  flight  conditions.  A  further  reduction  in  the  model  wing  cross 
sectional  properties  could  possibly  result  in  matching  the  wing  twist  through  the  varia¬ 
tion  in  load  factor. 


Figure  62  also  indicates  that  the  airplane  wing  twists  0.  38  degree  more  than  that  of 
the  model  at  a  load  factor  of  2.5.  The  additional  twist  means  that  the  outboard  airfoil 
sections  are  at  a  lower  angle  of  attack  than  the  corresponding  sections  on  the  wind 
tunnel  model.  The  increased  twist  results  in  a  reduced  lift,  drag,  and  wing  pitching 
moment  relative  to  the  wind  tunnel  model.  The  resulting  longitudinal  incremental  aero¬ 
dynamic  coefficients  from  the  differential  tip  twist  are  presented  in  Figure  63  as  a 
function  of  increased  load  factor. 


Figure  64  illustrates  the  aerodynamic  increments  as  a  function  of  the  model  angle  of 
attack.  The  jig  twisted  design  point  corresponds  to  the  point  on  Figure  64  at  6.0  degrees 
angle  of  attack.  For  this  condition,  the  airplane  is  flying  at  a  mid-cruise  weight  of 
234,000  pounds,  an  altitude  of  36,000  feet  at  Mach  0.  98,  and  a  Cl  of  0.  32.  The  incre¬ 
mental  drag  coefficient  varies  from  zero  at  the  cruise  condition  (nz=  1.  0)  to  0.  0012  at  a 
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Figure  62  .  Effect  of  Load  Factor  on  the  ATT 
Wing  Tip  Twist  at  Mach  0.  98 
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Figure  63.  Estimated  Effect  of  Load  Factor  on  the 
Incremental  Aerodynamic  Data  for  the 
ATT  Study  Configuration  at  Mach  0.  98 
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load  factor  of  2. 5  at  11. 4  degrees  angle  of  attack.  A  drag  increment  of  0.  0012  results 
in  a  drag  force  of  about  37  pounds  at  tunnel  conditions  simulating  full  scale  Reynolds 
number. 

The  model  has  been  designed  to  be  supported  in  the  HIRT  facility  on  a  sting  with  an 
internal  strain  gage  balance  used  to  measure  the  aerodynamic  loads.  The  balance 
designed  for  the  ATT  model  (Reference  5)  has  an  estimated  maximum  normal  force 
capacity  of  27,000  pounds  with  an  estimated  maximum  axial  force  capacity  of  2,700 
pounds.  The  accuracy  of  the  balance  has  been  estimated  at  about  0.  5  per  cent  of  the 
full-scale  load.  This  results  in  an  accuracy  of  about  ±13.  5  pounds  of  axial  force  with 
a  full-scale  loading.  Rotating  the  drag  force  increment  due  to  aeroelastic  differences 
with  the  accompanying  lift  force  increment  through  the  model  angle  of  attack  results 
in  an  axial  force  of  22  pounds.  This  estimated  drag  difference  is  within  the  accuracy 
of  the  balance;  however,  at  reduced  angles  of  attack,  an  incremental  drag  coefficient 
of  0.  0010  for  the  ATT  model  represents  the  balance  accuracy.  Below  this  level,  the 
aerodynamic  increments  cannot  be  measured  accurately,  thus  the  aeroelastic  defor¬ 
mations  will  not  significantly  affect  the  test  data  up  to  10.  5  degrees  angle  of  attack  for 
tunnel  conditions  at  full  scale  Reynolds  number  and  Mach  0.  98. 

5.3.2  Characteristics  with  Reduced  Wing  Cross  Sectional  Properties 

In  the  preceding  discussion,  the  model  wing  deflection  data  indicated  that  a  further 
reduction  in' the  model  wing  cross-sectional  properties  could  result  in  matching  the 
wing  twist  through  the  variation  in  load  factor.  A  model  condition  has  been  analyzed 
with  the  cross-sectional  properties  arbitrarily  reduced  from  the  65  percent  section 
to  a  52  percent  section.  The  initial  65  percent  section  was  achieved  by  allowing  the 
aft  35  percent  of  the  wing  to  be  removable.  The  removable  section  would  allow  access 
to  possible  wing  pressure  tubing.  By  reducing  the  wing  section  properties  to  52  per¬ 
cent,  the  model  wing  will  become  more  flexible.  The  increased  flexibility  results  in 
the  model  wing  twisting  more  than  the  aircraft  wing  for  increased  load  factors,  as 
shown  in  Figure  65 .  The  airplane  wing  twists  1.  7  degrees  per  g,  while  the  model 
twists  1.  9  degrees  per  g  under  the  HIRT  test  conditions.  Figure  65  indicates  that 
the  model  wing  twists  0. 45  degree  more  than  the  airplane  at  a  load  factor  of  2.  5. 

The  additional  model  wing  twist  results  in  test  data  with  a  reduced  lift  and  drag,  and 
increased  nose  up  pitching  moment  relative  to  the  airplane.  The  resulting  longitudinal 
aerodynamic  increments  from  the  differential  tip  twist  are  shown  in  Figure  66  as  a 
function  of  load  factor.  Figure  67  illustrates  the  estimated  difference  in  drag  due  to 
lift  as  a  function  of  model  angle  of  attack.  For  a  load  factor  of  2.5,  which  corresponds 
to  a  model  angle  of  attack  of  11.3  degrees,  the  incremental  drag  coefficient  was  0.0016. 
For  this  condition,  the  drag  increments,  lift  increments,  and  pitching  moment  incre¬ 
ments  must  be  added  to  the  model  test  data. 
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Figure  65.  Effect  of  Load  Factor  on  the  ATT  Wing  Tip  Twist 
at  Mach  0.  98  (52%  Model  Wing  Section) 
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Figure  66.  Estimated  Effect  of  Load  Factor  on  the  Incremental  Aerodynamic  Data 
for  the  ATT  Study  Configuration  at  Mach  0.  98  (52%  Model  Wing  Section) 
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Figure  67 .  Estimated  Effect  of  Angle  of  Attack  on  the  Incremental 
Aerodynamic  Data  for  the  ATT  Study  Configuration 
at  Mach  0.  98  (52%  Model  Wing  Section) 
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The  data  presented  for  the  two  cross-sectional  wing  properties  indicate  that,  for  this 
transport-type  airplane,  a  possible  cross-sectional  property  of  the  model  wing  could 
be  determined  and  designed  into  the  model,  which  would  cause  the  model  deformation 
to  be  the  same  as  that  of  the  flight  vehicle.  The  stress  analysis  for  the  reduced  cross- 
sectional  properties  of  the  model  wing  was  checked  and  found  to  be  strong  enough  for 
a  portion  of  the  flight  envelope  but  not  for  all  critical  points  within  the  envelope. 

5.3.3  Characteristics  at  Mach  0.80  and  25,  000  Feet 

The  previous  discussion  has  considered  the  ATT  aircraft  operating  at  a  high  transonic 
Mach  number  and  at  a  high  altitude.  Four  additional  points  have  been  analyzed  within 
the  flight  envelope,  based  on  the  model  wing  jig  twist  being  held  the  same  as  the  high 
altitude,  high  transonic  Mach  number  design  condition. 

For  a  flight  condition  of  25,  000  feet  and  Mach  0.80,  the  aerodynamic  analysis  resulted 
in  the  wing  twist  data  for  the  model  and  airplane  as  shown  in  Figure  68  for  increasing 
load  factors.  The  data  in  Figure  68  indicate  that  the  aircraft  wing  twists  2.0  degrees 
per  g,  while  the  model  twists  1. 75  degrees  per  g  under  the  HIRT  test  conditions.  At 
a  load  factor  of  2.5,  the  incremental  wing  twist  was  0.5  degree,  with  the  airplane 
twisting  more  than  the  model  wing.  Longitudinal  aerodynamic  increments  are  presented 
in  Figure  69  as  a  function  of  the  load  factor  and  as  a  function  of  the  model  angle  of 
attack  are  in  Figure  70.  At  a  load  factor  of  2.  5,  the  differential  drag  coefficient  was 
0.0012,  which  results  in  a  drag  force  of  about  39. 5  pounds  at  tunnel  conditions  simula¬ 
ting  full  scale  Reynolds  number.  The  balance  accuracy  was  13.5  pounds  of  axial  force. 
Rotating  the  drag  increment  through  the  angle  of  attack  results  in  an  axial  force  of  51 
pounds,  which  is  within  the  accuracy  of  the  balance.  However,  with  reduced  angle  of 
attack,  an  incremental  drag  coefficient  of  0.0007  represents  the  balance  accuracy  for 
the  tunnel  conditions  in  this  case. 

5.3.4  Characteristics  of  Mach  0,95  and  20,000  Feet 

For  a  flight  condition  of  20,000  feet  and  Mach  0.95,  the  aeroelastic  analysis  resulted 
in  the  wing  twist  data  shown  in  Figure  71  for  increasing  load  factor.  The  data  indicate 
that  the  aircraft  wing  twists  1. 75  degrees  per  g,  with  the  model  wing  twisting  1.40 
degrees  per  g  under  the  HIRT  test  conditions.  In  level  flight  at  a  of  0. 19,  the  air¬ 
plane  wing  tip  twists  0.60  degree  more  than  the  model  wing.  At  a  load  factor  of  2.5, 
this  airplane  wing  twists  1. 10  degrees  more  than  the  model  wing.  The  resulting  aero¬ 
dynamic  increments  as  a  function  of  the  load  factor  are  presented  in  Figure  72,  with 
the  data  presented  in  Figure  73  as  a  function  of  the  model  angle  of  attack.  At  a  load 
factor  of  2.5,  the  incremental  drag  coefficient  was  0.0015,  which  results  in  a  drag 
force  of  about  74  pounds  at  tunnel  conditions  simulating  full  scale  Reynolds  number. 
Rotating  the  drag  increment  through  the  angle  of  attack  results  in  an  axial  force  of 
22  pounds.  This  is  again  within  the  accuracy  of  the  balance;  however,  with  reduced 
angles  of  attack  an  incremental  drag  coefficient  of  0.0014  represents  the  balance 
accuracy  for  the  tunnel  conditions  in  this  case. 
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Figure  68  .  Effect  of  Load  Factor  on  the  ATT 
Wing  Tip  Twist  at  Mach  0.  8 
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Figure  69  .  Estimated  Effect  of  Load  Factor  on  the  Incremental  Aero¬ 
dynamic  Data  for  the  ATT  Study  Configuration  at  Mach  0.  8 
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Figure  71 .  Effect  of  Load  Factor  on  the  ATT 
Wing  Tip  Twist  at  Mach  0.  95 
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Figure  72  .  Estimated  Effect  of  Load  Factor  on  the  Incremental  Aerodynamic 
Data  for  the  ATT  Study  Configuration  at  Mach  0.  95 
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Figure  73  .  Estimated  Effect  of  Angle  of  Attack  on  the  Incremental  Aerodynamic 
Data  for  the  ATT  Study  Configuration  at  Mach  0.  95 


5.3.5  Characteristics  at  Mach  0.  60  and  10,  000  Feet 

A  flight  condition  at  10, 000  feet  and  Mach  0.  6  was  analyzed.  The  wing  tip  twist  data 
are  presented  in  Figure  74  .  The  data  indicate  that  the  aircraft  wing  twists  2.  0  degrees 
per  g,  with  the  model  wing  twisting  1.  70  degrees  per  g  under  the  HIRT  test  conditions. 
At  a  load  factor  of  2.  5,  the  incremental  tip  twist  was  0.  7  degree,  with  the  aircraft  wing 
deflecting  more  than  the  model  wing.  Aerodynamic  increments  are  presented  in  Fig¬ 
ure  75  as  a  function  of  load  factor  and  in  Figure  76  as  a  function  of  the  model  angle  of 
attack.  For  this  condition,  this  drag  increment  varies  from  0.  0003  at  the  level  flight 
condition  at  a  CL  of  0. 32,  to  0.  0015  at  a  load  factor  of  2.  5  at  a  Cl  of  0.  80.  The  incre¬ 
mental  drag  coefficient  of  0.  0015  results  in  a  drag  force  of  48  pounds  at  the  tunnel  con¬ 
ditions  simulating  full  scale  Reynold  number  and  at  Mach  0.  6.  Rotating  the  drag  in¬ 
crement  through  the  angle  of  attack  results  in  an  axial  force  of  -55  pounds.  This  is 
again  within  the  accuracy  of  the  balance;  however,  with  a  reduced  angle  of  attack,  an 
incremental  drag  coefficient  of  0.  0009  represents  the  balance  accuracy  for  the  tunnel 
conditions  in  this  case. 


Figure  77  illustrates  selected  data  of  the  model  wing  tip  twist  as  a  function  of  Mach 
number  at  a  constant  dynamic  pressure.  The  data  indicate  only  a  small  effect  of  Mach 
number  on  the  wing  tip  twist. 
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Figure  74  .  Effect  of  Load  Factor  on  the  ATT 
Wing  Tip  Twist  at  Mach  0.  60 
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Figure  75  .  Effect  of  Load  Factor  on  the  Incremental  Aerodynamic 
Data  for  the  ATT  Study  Configuration  at  Mach  0.  60 
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Figure  76.  Estimated  Effect  of  Angle  of  Attack  on  the  Incremental  Aero¬ 
dynamic  Data  for  the  ATT  Study  Configuration  at  Mach  0. 60 


Figure  77.  Effect  of  Mach  Number  on  the  ATT  Wing  Tip  Twist 
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5.3.6  Summary 

The  ATT  transport  model  used  in  this  study  was  considered  to  be  designed  as  an  all 
steel  model  with  35  percent  of  the  wing  section  not  considered  as  contributing  to  the 
bending  and  torsional  stiffness  characteristics  because  of  the  inclusion  of  wing- 
pressure  tube  routing  channels.  The  aeroelastic  analysis  of  the  wing  resulted  in  less 
model  wing  twist  and  vertical  deflection  than  the  airplane  wing  at  equivalent  Reynolds 
number  conditions  as  shown  in  Figures  68  and  56.  In  the  aerodynamic  analysis, 
vertical  deflection  or  the  effect  of  a  nonplanar  wing  was  not  considered  directly  but 
was  considered  insofar  as  including  the  primary  effect  of  vertical  deflection  on  twist. 
Further  effect  of  a  nonplanar  wing  on  the  aerodynamic  characteristics  was  considered 
to  be  very  small. 

The  deformation  differences  between  the  model  wing  and  the  airplane  do  not  result  in 
significant  differences  in  the  estimated  model  and  airplane  aerodynamics.  For  in¬ 
stance,  at  the  design  cruise  condition  (M=0.98,  Alt.  =36,000  ft)  and  at  the  maximum 
load  factor  (2.5)  the  drag  of  the  model  was  found  to  be  0.0012  greater  than  the  air¬ 
plane.  For  this  test  condition,  the  model  force  balance  accuracy  in  drag  was  considered 
to  be  about  ±0.0010.  Thus,  the  balance  accuracy  is  of  as  much  concern  as  the  aero¬ 
elastic  effects  on  the  aerodynamic  data  for  this  configuration  since  at  least  approximate 
corrections  can  be  made  for  the  aerodynamic  differences  due  to  aeroelasticfty  if  the 
deflection  are  lmown. 
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5.  4  LONGITUDINAL  AEROELASTIC  CHARACTERISTICS  OF  THE  F-lll 

Aeroelastic  deformations  of  the  F-lll  wing  were  analyzed  for  several  flight  conditions 
of  the  HIRT  model  through  the  use  of  the  modified  Weissinger  L-method.  For  the  F-lll 
model,  the  deflected  wing  shape  of  the  aircraft  in  level  flight  (n  =  1)  at  10,000  feet  and 
Mach  0.  90  was  used  to  determine  the  model  wing  jig  twist.  The  resulting  jig  twist,  as 
shown  in  Figure  78  ,  resulted  in  matching  the  model  wing  deflection  with  the  aircraft 
wing  deflection.  The  model  wing  design  (Reference  3)  was  assumed  to  be  constructed 
of  solid  steel  without  cutouts  in  the  wing  cross-sectional  characteristics. 

5.4.1  Characteristics  at  Mach  0.90  and  10,  000  Feet 

Figure  78  illustrates  a  comparison  of  the  airplane  wing  twist  with  that  of  the  model  for 
HIRT  test  conditions  simulating  the  full  scale  Reynolds  number  at  10, 000  feet  and  Mach 
0.90.  Figure  79  illustrates  the  comparison  of  the  airplane  wing  tip  twist  with  that  of 
the  model  for  increasing  load  factors.  The  data  indicate  that  the  model  wing  twists  0.7 
degree  per  g  under  the  HIRT  test  conditions,  with  the  airplane  wing  twisting  0.6  degree 
per  g  voider  flight  conditions.  The  model  wing  twists  about  0.5  degree  more  than  the 
airplane  wing  at  a  load  factor  of  5.0.  For  this  case,  the  model  is  more  flexible  under 


Figure  78  .  Comparison  of  the  F-lll  Model  Wing  Twist  with  the 
F-lll  Airplane  Wing  for  Level  Flight  at  Mach  0.  90 
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Figure  79.  Effect  of  Load  Factor  on  the  F-lll  Wing- 
Tip  Twist  at  Mach  0.  90  and  Full-Scale 
Reynolds  Number  (Altitude  =  10,000  ft) 

the  HIRT  test  conditions  than  the  airplane  under  flight  conditions.  From  this  condition, 
one  of  the  possible  methods  of  improving  the  resulting  model  aeroelastic  deforma¬ 
tions  would  be  the  tunnel  dynamic  pressure.  This  condition  is  discussed  later  in  this 
section. 

The  additional  twist  experienced  by  the  model  wing  means  that  the  outboard  airfoil 
sections  are  at  a  lower  angle  of  attack,  which  results  in  a  decrease  in  lift  and  drag 
and  more  nose-up  pitching  moment  relative  to  the  airplane. 

The  longitudinal  aerodynamic  increments  resulting  from  the  additional  tip  twist  are 
presented  in  Figure  80  as  a  function  of  the  aircraft  load  factor.  Figure  81  illustrates 
the  aerodynamic  increments  as  a  function  of  the  model  angle  of  attack. 

The  jig  twist  design  point  corresponds  to  the  point  on  Figure  81  at  2.4  degrees  angle 
of  attack.  For  this  condition,  the  airplane  is  flying  at  a  cruise  weight  of  70,000  pounds, 
an  altitude  of  10, 000  feet,  a  speed  of  Mach  0.  90,  and  a  Cj,  of  0.  16.  The  incremental 
drag  coefficient  varies  from  zero  at  the  matched  condition  to  0.  0010  at  a  load  factor  of 
5.  0  at  10. 3  degrees  angle  of  attack.  A  drag  increment  of  0.  0010  results  in  a  drag 
force  of  about  29.4  pounds  at  tunnel  conditions  simulating  full  scale  Reynolds  number. 

The  model  has  been  designed  to  be  supported  in  the  HIRT  facility  on  a  sting  with  an 
internal  strain  gage  balance  used  to  measure  the  aerodynamic  loads.  The  balance 
design  for  the  F-lll  model  (Reference  5)  has  an  estimated  maximum  normal  force 
capacity  of  37,000  pounds,  with  an  estimated  maximum  axial  force  of  3,700  pounds. 
Balance  accuracy  has  been  estimated  to  be  about  0.  5  percent  of  the  full  scale  load. 

This  results  in  an  accuracy  of  about  18.  5  pounds  of  axial  force  with  a  maximum  axial 
force  loading. 

Rotating  the  drag  force  increment  due  to  aeroelastic  differences  with  the  accompanying 
lift  force  through  the  angle  of  attack  results  in  an  axial  force  of  23.  9  pounds.  This  is 
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Figure  80 .  Estimated  Effect  of  Load  Factor  on  the 
Incremental  Aerodynamic  Data  for  the 
F-lll  Study  Configuration  at  Mach  0. 90 

within  the  accuracy  of  the  balance;  however,  with  reduced  angles  of  attack,  an  incre¬ 
mental  drag  coefficient  of  0.  0006  represents  the  balance  accuracy.  Below  this  level, 
the  aerodynamic  drag  increment  due  to  aeroelastic  effects  cannot  be  measured  accurately, 
thus  the  test  data  will  not  be  affected  up  to  about  7.  0  degrees  angle  of  attack. 

5.  4.  2  Characteristics  at  75  Percent  Full  Scale  Test  Reynolds  Number 

For  certain  conditions  it  may  be  required  to  minimize  the  aerodynamic  increments  as 
presented  in  Figure  80.  To  achieve  this ,  the  HIRT  facility  offers  unique  methods  of 
dynamic  pressure  reduction.  First  would  be  simply  a  reduction  in  the  test  Reynolds 
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number;  secondly  would  be  a  decrease  in  the  tunnel  operating  temperature,  which 
would  allow  full  scale  Reynolds  number  testing  at  a  lower  dynamic  pressure.  Three 
test  conditions  have  been  considered  and  analyzed  for  the  study  design  point  on  the 
F-lll  HIRT  model. 

For  a  25-percent  reduction  in  the  test  Reynolds  number  (Retest  =  0*75  Reairpiane)  > 

the  incremental  wing  twist  was  reduced  from  0. 5  to  0. 25  degree  at  a  load  factor  of 
5.  0.  Figure  82  illustrates  the  wing  tip  twist  comparison  for  a  25-percent  reduction 
in  the  full  scale  test  Reynolds  number  as  a  function  of  the  load  factor.  For  this  con¬ 
dition,  the  model  wing  does  not  twist  as  much  as  the  airplane  wing  under  flight  condi¬ 
tions.  The  outboard  sections  of  the  model  wing  are  at  a  slightly  greater  angle  of 
attack  than  the  airplane  wing.  This  results  in  the  aerodynamic  increments  as  shown 
in  Figure  83.  The  test  data  has  been  estimated  to  indicate  a  higher  drag  and  lift  and 
greater  nose-up  pitching  moment  than  the  airplane  would  indicate.  Figure  84  shows 
the  aerodynamic  increments  as  a  function  of  the  model  angle  of  attack.  For  a  25- 
percent  reduction  in  the  full  scale  test  Reynolds  number,  the  results  indicate  that  the 
aerodynamic  corrections  are  smaller  than  for  tests  performed  at  full  scale  test  Reynolds 
number. 

5.  4.  3  Characteristics  at  50  Percent  Full  Scale  Test  Reynolds  Number 

For  a  50  percent  reduction  in  the  test  Reynolds  number,  the  estimated  difference  in 
the  model  and  the  airplane  aerodynamic  characteristics  became  significantly  greater. 

The  tip  twist  is  presented  in  Figure  85  as  a  function  of  load  factor  for  the  50  percent 
reduced  test  Reynolds  number  case,  with  the  airplane  at  flight  conditions.  The  differ¬ 
ential  wing  tip  twist  angle  becomes  1.  45  degrees  at  a  load  factor  of  5.  0.  The  incre¬ 
mental  aerodynamic  terms ,  as  presented  in  Figures  86  and  87 ,  become  greater  in 
magnitude  than  those  in  Figures  80  and  81  at  full  scale  test  Reynolds  number. 

Thus,  fonthe  F-lll  study  configuration,  a  small  reduction  in  test  Reynolds  number 
will  improve  the  model  simulation  in  terms  of  lift  and  drag  due  to  lift  and  pitching 
moment. 

5.4.4  Characteristics  at  Full  Scale  Test  Reynolds  Number  with  Reduced  Test  Section 
Temperature 


The  preceding  discussions  on  reducing  the  test  Reynolds  number  to  help  the  model 
loading  conditions  is  one  method  of  minimizing  the  aeroelastic  effects  on  the  test  data. 

A  second  method  proposed  for  the  HERT  facility  is  to  control  the  tunnel  air  temperature. 
This  allows  a  reduction  in  the  tunnel  air  temperature  for  testing  at  full  scale  Reynolds 
number  with  a  reduced  dynamic  pressure.  For  this  tunnel  condition,  a  data  point  was 
generated  on  the  F-lll  study  configuration  based  on  a  temperature  drop  from  300° K  to 
240°  K.  The  tunnel  characteristics  from  reference  13  were  used  to  determine  the  tunnel 
dynamic  pressure  at  the  240° K  air  temperature.  The  results  produced  a  dynamic  pres¬ 
sure  identical  to  that  for  the  25-percent  reduction  in  the  test  Reynolds  number  case  dis¬ 
cussed  in  Section  5.  4. 2. 
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Figure  82.  Effect  of  Load  Factor  on  the  F-lll  Wing  Tip  Twist 
at  Mach  0.  90  and  75%  Full  Scale  Reynolds  Number 
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Figure  83.  Estimated  Effect  of  Load  Factor  on  the  Incremental 
Aerodynamic  Data  for  the  F-lll  Study  Configuration 
at  Mach  0.  90  and  75%  Full  Scale  Reynolds  Number 
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Figure  84.  Estimated  Effect  of  Angle  of  Attack  on  the  Incremental 
Aerodynamic  Data  for  the  F-lll  Study  Configuration 
at  Mach  0.  90  and  75%  Full  Scale  Reynolds  Number 
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Figure  85.  Effect  of  Load  Factor  on  the  F-lll  Wing  Tip  Twist 
at  Mach  0.  90  and  50%  Full  Scale  Reynolds  Number 
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Figure  86.  Estimated  Effect  of  Load  Factor  on  the  Incremental 
Aerodynamic  Data  for  the  F-lll  Study  Configuration 
at  Mach  0.  90  and  50%  Full  Scale  Reynolds  Number 
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Figure  87.  Estimated  Effect  of  Angle  of  Attack  on  the  Incremental 
Aerodynamic  Data  for  the  F-lll  Study  Configuration 
at  Mach  0.  90  and  50%  Full  Scale  Reynolds  Number 
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The  tip  twist  for  the  240°  K  tunnel  temperature  condition  is  presented  in  Figure  88. 

The  incremental  aerodynamic  terms  are  presented  in  Figures  89  and  90.  For  this 
case,  the  reduction  from  ambient  conditions  at  300°  K  to  the  cold  condition  at  240°  K 
results  in  aerodynamic  corrections  that  are  smaller  than  for  tests  performed  at  the 
300°  K.  Thus,  for  the  F-lll  study  configuration,  a  reduction  in  the  tunnel  temperature 
would  result  in  improved  model  simulation  in  terms  of  the  aerodynamic  increments  in 
lift,  drag,  and  pitching  moment  without  a  reduction  in  the  test  Reynolds  number. 

The  elastic  twist  data  presented  in  this  section,  together  with  the  cases  at  reduced 
test  Reynolds  number,  indicate  an  incremental  tip  twist  variation  from  0.  5  to  1.  45 
degrees  at  a  load  factor  of  5.  0.  The  proposed  wing  twist  measuring  system  for  the 
HIRT  facility  will  definitely  be  of  great  importance  to  the  interpretation  of  the  test 
data  over  the  ranges  of  variables  discussed  in  Sections  5.  4. 1  through  5.  4. 3.  The  data 
indicate  the  sensitivity  to  the  tunnel  dynamic  pressure  due  to  a  reduction  in  the  test 
Reynolds  number  or  the  tunnel  operating  temperature. 

5.  4.  5  Characteristics  at  Mach  0.  90  and  30,  000  Feet 

A  case  at  Mach  0.  90  and  30, 000  feet  was  analyzed  for  the  F-lll  study  configuration. 

The  tip  twist  is  compared  for  the  HIRT  model  and  the  airplane  for  flight  load  factors 
in  Figure  91.  The  data  indicate  that  the  model  wing  twists  0.  9  degree  per  g  under  the 
HIRT  test  conditions.  The  airplane  wing  twists  0.  7  degree  per  g  under  flight  conditions. 
The  model  wing  twists  0.  4  degree  more  than  the  airplane  wing  at  a  load  factor  of  2.  0. 
The  incremental  aerodynamic  terms  are  presented  in  Figure  92  as  a  function  of  load 
factor  and  in  Figure  93  as  a  function  of  the  model  angle  of  attack. 

5. 4.  6  Characteristics  at  Mach  0.90  and  20,000  Feet 

A  flight  condition  at  Mach  0,  90  and  20,000  feet  was  analyzed  for  the  F-lll  airplane  and 
the  HIRT  model  under  the  HIRT  test  conditions  corresponding  to  full  scale  Reynolds 
number  and  Machnumber.  The  tip  twist  is  compared  for  the  model  and  the  airplane  in 
Figure  94.  The  data  indicate  that  the  model  wing  twists  0.  80  degree  per  g  under  the 
HIRT  test  conditions,  while  the  airplane  wing  twists  0.  53  degree  per  g  under  the  flight 
conditions.  The  model  wing  twists  0.  65  degree  more  than  the  airplane  wing  at  a  load 
factor  of  4.  0.  The  incremental  aerodynamic  terms  that  result  from  the  differential 
tip  twist  are  presented  in  Figure  95  as  a  function  of  the  load  factor.  Figure  96  illus¬ 
trates  the  incremental  terms  as  a  function  of  the  model  angle  of  attack. 

5.  4.  7  Characteristics  at  Mach  0.  70  and  10,000  Feet 

A  flight  condition  at  Mach  0.  70  and  10,000  feet  was  analyzed  for  the  F-lll  airplane 
and  the  model  under  the  HIRT  test  conditions  corresponding  to  full  scale  Reynolds 
number  and  Mach  number.  The  tip  twist  is  compared  for  the  model  and  the  airplane 
in  Figure  97.  The  data  indicate  that  the  model  wing  twists  0.  80  degree  per  g  under 
the  HIRT  test  conditions  while  the  airplane  wing  twists  0.  50  degree  per  g  under  flight 
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Figure  88.  Effect  of  Load  Factor  on  the  F-lll  Wing  Tip  Twist 
at  Mach  0.  90  and  Full  Scale  Reynolds  Number  with 
the  HTRT  Tunnel  Temperature  Reduced  to  240°K 


LOAD  FACTOR  (g) 

Figure  89.  Estimated  Effect  of  Load  Factor  on  the  Incremental 
Aerodynamic  Data  for  the  F-lll  Study  Configuration 
at  Mach  0.  90  and  Full  Scale  Reynolds  Number  with 
the  HIRT  Tunnel  Temperature  Reduced  to  240°K 
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Figure  90.  Estimated  Effect  of  Angle  of  Attack  on  the  Incremental 
Aerodynamic  Data  for  the  F-lll  Study  Configuration 
at  Mach  0.  90  and  Full  Scale  Reynolds  Number  with 
the  HIRT  Tunnel  Temperature  Reduced  to  240°K 


Figure  91.  Effect  of  Load  Factor  on  the  F-lll  Wing 
Tip  Twist  at  Mach  0.  90  and  Full  Scale 
Reynolds  Number  (Altitude  =  30,000  ft) 
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Figure  92.  Estimated  Effect  of  Load  Factor  on  the  Incremental 
Aerodynamic  Data  for  the  F-lll  Study  Configuration 
at  Mach  0.  90  and  Full  Scale  Reynolds  Number 
(Altitude  =  30,000  ft) 

conditions.  The  model  wing  twists  0.  60  degree  more  than  the  airplane  wing  at  a  load 
factor  of  3.  4.  The  incremental  aerodynamic  terms  that  result  from  the  differential  tip 
twist  are  presented  in  Figure  98  as  a  function  of  the  load  factor.  Figure  99  illustrates 
the  incremental  terms  as  a  function  of  the  model  angle  of  attack. 

5.  4.  8  Tip  Twist  Characteristics  as  a  Function  of  Dynamic  Pressure 
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Figure  100  illustrates  the  wing  tip  twist  as  a  function  of  dynamic  pressure  for  the  F-lll 
HIRT  model  at  a  constant  Mach  number.  For  the  condition  of  level  flight,  the  wing  tip 
twist  remains  relatively  constant  with  increasing  dynamic  pressure.  At  load  factors 
greater  than  1.  0,  the  tip  twist  decreases  with  increasing  dynamic  pressure. 
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Figure  93.  Estimated  Effect  of  Angle  of  Attack  on  the  Incremental  Aero¬ 
dynamic  Data  for  the  F-lll  Study  Configuration  at  Mach  0.  90 
and  Full  Scale  Reynolds  Number  (Altitude  =  30,000  ft) 
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Figure  94.  Effect  of  Load  Factor  on  the  F-111  Wing 
Tip  Twist  at  Mach  0.  90  and  Full  Scale 
Reynolds  Number  (Altitude  =  20,000  ft) 
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Figure  95.  Estimated  Effect  of  Load  Factor  on  the  Incremental  Aero¬ 
dynamic  Data  for  the  F-111  Study  Configuration  at  Mach  0.  90 
and  Full  Scale  Reynolds  Number  (Altitude  =  20,000  ft) 
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Figure  98.  Estimated  Effect  of  Load  Factor  on  the  Incremental  Aerody¬ 
namic  Data  for  the  F-lll  Study  Configuration  at  Mach  0.  70 
and  Full  Scale  Reynolds  Number  (Altitude  =  10,000  ft) 


97 


AEDC-TR-75-62 


0.0015. 
0.0010 
ACd  0.0005 
0 

-0.0005 

“0.0040 

-0.0030 

-0.0020 

AC 

m 

-0.0010 

0 

0.0010 
0.0150 
0.0100 
acl  0.0050 
0 

-0.0050 


n  =  3.0 

2  I 


n  =  2 . 0 

^  I 


WING  SWEEP  =  50° 

—  F-  111  AIRPLANE:  M-O. 70.  q=498  psf,  ALT-1 0, 000  ft  - 
F-111  HI  RT  MODEL,  M-0.70,  q=5,256psf,  TT=300°K 
FULL  SCALE  REYNOLDS  NUMBER 


a  (deg) 


Figure  99.  Estimated  Effect  of  Angle  of  Attack  on  the  Incremental  Aero¬ 
dynamic  Data  for  the  F-111  Study  Configuration  at  Mach  0.  70 
and  Full  Scale  Reynolds  Number  (Altitude  =  10,000  ft) 
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Figure  100.  Effect  of  Tunnel  Dynamic  Pressure  on  the  F-lll  HURT  Model 
Wing  Tip  Twist  at  Mach  0.90  and  Full  Scale  Reynolds  Number. 
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5.4.9  Summary 

The  F-lll  model  used  in  this  study  was  considered  to  be  designed  as  an  all  steel  model 
with  solid  steel  wings .  The  aeroelastic  analysis  of  the  wing  resulted  in  greater  model 
wing  twist  and  vertical  deflection  than  the  airplane  wing  at  equivalent  Reynolds  number 
conditions  as  shown  in  Figures  79  and  40.  In  the  aerodynamic  analysis,  vertical  de¬ 
flection  of  the  wing  was  not  considered  directly  but  was  considered  insofar  as  including 
the  primary  effect  of  vertical  deflection  on  twist.  The  further  effect  of  a  nonplanar 
wing  on  aerodynamic  characteristics  was  considered  to  be  very  small.  The  deforma¬ 
tion  differences  between  the  model  wing  and  the  airplane  do  not  result  in  significant 
differences  in  the  estimated  model  and  airplane  aerodynamics.  For  instance,  at  the 
design  cruise  condition  (M=0.90,  Alt.  =10,000  ft)  and  at  a  high  load  factor  (5.0)  the 
drag  of  the  model  was  0.0010  less  than  the  airplane.  For  this  test  condition,  the 
model  force  balance  accuracy  in  drag  was  considered  to  be  about  ±0.0006.  Thus,  as 
was  stated  for  the  ATT  model,  the  balance  accuracy  is  of  as  much  concern  as  the 
aeroelastic  effects  on  the  aerodynamics  data. 

The  effect  of  mismatching  test  Reynolds  number  and  reduced  tunnel  temperatures 
were  analyzed  for  the  F-lll  model.  The  data  at  a  level  of  75%  full  scale  Reynolds 
number  as  shown  in  Figures  82,  83  and  84  decreased  the  model/airplane  deformation 
differences  and  also  the  incremental  aerodynamic  terms.  For  a  50%  full  scale 
Reynolds  number  test  condition,  the  data  as  shown  in  Figure  85,  86  and  87  indicates 
increased  model/airplane  deformation  differences  and  increased  incremental  aero¬ 
dynamic  terms.  This,  only  slight  mismatching  of  Reynolds  number  will  bring  the 
model/airplane  deformations  into  agreement  along  with  the  aerodynamic  terms. 

The  effect  of  reducing  tunnel  storage  temperature  from  300°  K  to  240°  K  was  analyzed 
for  the  F-lll  model.  The  data  indicates  smaller  model/airplane  deformation  and 
aerodynamic  differences  can  be  achieved  through  the  tunnel  storage  temperature 
control.  For  this  configuration,  where  the  model  deformations  are  greater  than  the 
airplane  deformations  at  equivalent  Reynolds  number  conditions,  the  control  of  tunnel 
temperature  may  be  the  only  method  of  improving  the  model/airplane  deformations  and 
thus  the  aerodynamic  terms. 
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SECTION  VI 
CONCLUSIONS 

This  study  indicates  that  for  the  two  study  configurations,  a  model  wing  jig  twist  can 
be  selected  such  that  the  elastic  deformation  of  the  model  wing  will  be  similar  to  the 
airplane  for  a  wide  range  of  conditions  within  the  flight  envelope.  The  resulting  aero- 
elastic  differences  (i.e.  twist)  between  the  model  and  the  airplane  at  an  equivalent 
Reynolds  number  will  not  result  in  significant  effects  on  the  aerodynamic  test  data. 


The  HIRT  facility  will  provide  more  accurate  test  data  because  of  the  model  flexibility 
than  current  "rigid"  model  test  techniques  will  provide.  The  proposed  model  defor¬ 
mation  measuring  system  will  be  essential  since  accurate  measurement  of  the  model 
deflection  will  provide  information  for  correct  interpretation  of  the  test  data. 


The  study  further  indicated  that  three  options  are  available  to  produce  close  similarity 
of  the  model/airplane  wing  deformation  over  a  broad  range  of  operating  conditions. 
These  include  mismatching  Reynolds  number,  variation  in  the  tunnel  operating 
temperature,  and  construction  of  more  than  one  wing.  For  configurations  where  model 
aeroelastic  deformations  are  excessive,  dynamic  pressure  reduction  through  a  reduc¬ 
tion  in  the  test  Reynolds  number  reduces  the  model  wing  deformations  so  that  geometric 
similarity  between  the  model  and  airplane  may  be  attained.  Reduction  in  the  tunnel 
temperature  such  as  from  300°  K  to  240°  K,  reduces  the  model  wing  deformations  so 
that  geometric  model/airplane  similarity  is  attained  at  equivalent  Reynolds  numbers. 
For  cases  where  model  aeroelastic  deformations  are  excessive  after  reduction  in 
tunnel  dynamic  pressure  conditions,  the  fabrication  of  an  additional  wing  may  be  re¬ 
quired.  The  additional  wing  with  a  new  jig  twist  distribution  would  allow  improvement 
in  the  model/airplane  similitude  for  these  conditions. 


Finally,  based  on  the  HIRT  balance  study,  (Reference  5),  the  balance  accuracy  will 
be  as  significant  as  the  estimated  aerodynamic  differences  due  to  aeroelastic  differences 
between  a  model  and  the  airplane. 
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APPENDIX  A 

A  STUDY  OF  HIRT  MODEL  AEROELASTIC  CHARACTERISTICS  IN  REFERENCE  TO 
THE  AEROELASTIC  NATURE  OF  AN  AIRPLANE  COMPOSITE  WING 

The  main  part  of  this  study  was  concerned  with  determining  and  assessing  the  differ¬ 
ences  in  model/airplane  wing  deformations  and  the  effects  of  these  differences  on  the 
longitudinal  aerodynamic  characteristics.  Both  study  configurations  were  considered 
to  have  conventional  aluminum  structures  for  the  flight  article,  while  the  wind  tunnel 
models  were  constructed  from  steel. 

Future  aircraft  designers  are  likely  to  depart  from  total  aluminum  construction  and 
make  use  of  composite  materials  such  as  graphite /epoxy  or  boron/aluminum.  In  this 
appendix,  the  representation  of  a  full-scale  composite  wing  by  a  steel  wind  tunnel  model 
wing  is  analyzed  to  determine  the  wing  deformations  under  load  and  assess  the  effect  of 
these  differences  on  the  aerodynamics.  The  analysis  has  been  performed  for  the  F-lll 
configuration  with  the  wing  designed  to  an  expected  value  of  structural  properties  for 
an  anisotropic  (composite)  structure.  Section  A1  of  this  appendix  describes  the  aero- 
elastic  characteristics  of  the  selected  composite  wing.  Section  A2  describes  the  result¬ 
ing  incremental  aerodynamic  characteristics  of  the  study  configuration.  Section  A3 
describes  the  analytical  results  for  simulating  a  composite  wing  with  a  solid  steel  model 
wing.  Section  A4  states  conclusions  of  this  study. 

Al.  AEROELASTIC  CHARACTERISTICS  OF  A  COMPOSITE  WING  FOR  THE  F-lll 
AIRPLANE 

Aeroelastic  deformations  of  the  F-lll  wing  were  analyzed  for  one  flight  condition  using 
the  same  approach  and  methods  previously  identified  in  this  report.  The  wing  was 
analyzed  with  assumed  values  of  structural  properties,  El  and  GJ,  for  an  anisotropic 
structure,  Anisotropicity  causes  the  wing  to  be  more  flexible  in  accommodating  all  of 
the  applied  loads  but  leads  to  an  increase  in  the  complexity  of  the  structural  wing  design 
due  to  the  large  combinations  of  ply  orientation.  Ply  orientation  can  be  selected  to  pro¬ 
duce  positive  or  negative  elastic  wing  twist  in  contrast  with  a  wing  constructed  from 
isotropic  materials.  For  this  analysis,  three  locations  of  the  elastic  axis  have  been 
selected  and  examined. 

Table  A-l  depicts  the  structural  properties,  E  and  G,  for  steel,  aluminum,  boron/ 
aluminum  and  graphite/epoxy  materials  with  crossply  percentages  ranging  from  20  to 
100  percent  for  the  graphite /epoxy. 

Figure  A-l  illustrates  the  E/G  ratio  graphically  for  the  percentage  of  crossplies  in  a 
graphite/epoxy  structure.  Figure  A-2  illustrates  several  airplane  wing  stiffness  ratios 
that  have  been  fabricated  or  proposed  (Reference  1).  It  is  expected  that  future  wing 
designs  will  tend  toward  higher  E/G  ratios  since  anisotropic  structures  may  allow  a 
design  of  the  structure  to  obtain  the  best  compromise  of  the  flutter  and  structural 
divergence  characteristics  with  a  low  structural  weight.  The  advanced  fighter  attack 
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Table  A-l.  Material  Structural  Properties 


Material 

E  (millions  of  psi) 

G  (millions  of  psi) 

E/G 

Aluminum 

10.5 

4.0 

2.62 

Steel 

29.5 

11.0 

2.  63 

Boron/ Aluminum 

32.0 

8.0 

4.00 

Graphite/Epoxy 

20%  Crossplies 

14.0 

1.50 

9.33 

38%  Crossplies 

11.  8 

2.34 

5.05 

67%  Crossplies 

7.7 

2.90 

2.66 

80%  Crossplies 

5.5 

4.  00 

1.37 

100%  Crossplies 

1.9 

5.30 

0.36 
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Figure  A-l.  Characteristics  of  Graphite /Epoxy  Structure 
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SEMISPAN 

Figure  A-2„  Stiffness  Ratio  Characteristics  of  Existing 
Airplanes  and  Proposed  Designs 

wing  characteristics  (Reference  1)  illustrated  in  Figure  A-2  have  a  nominal  stiffness 
ratio  of  3,  0.  This  level  of  nominal  stiffness  was  selected  to  represent  the  composite 
wing  for  this  study.  The  stiffness  ratio  of  3.0  corresponds  to  a  graphite/epoxy  struc¬ 
ture  with  38  percent  of  the  laminates  crossplied  relative  to  a  0°  axis  located  along  the 
wing  center  spar.  The  crossplied  laminates  may  be  adjusted  to  any  orientation,  such 
as  ±45°  or  ±60°  relative  to  the  0°  axis.  By  appropriate  adjustment  of  the  proportions 
and  orientation,  virtually  any  design  requirement  can  be  met.  The  E/G  ratio  for  the 
graphite /epoxy  was  about  5.  0  and  the  i/j  ratio  for  the  airplane  wing  was  about  0.  6, 
which  resulted  in  the  stiffness  ratio  (El/GJ)  ratio  of  3.0. 

The  HIRT  wind  tunnel  model  wing  was  jig  twisted  and  analyzed  in  Section  IV  of  this 
report,  based  on  the  design  twist  and  isotropic  characteristics  of  the  F-lll  aluminum 
wing.  For  this  effort,  the  model  wing  twist  was  not  changed.  The  composite  wing  jig 
twist  distributions  for  three  locations  of  the  elastic  axis  was  determined  by  holding  the 
loaded  wing  twist  distribution  identical  to  the  aluminum  wing  for  level  flight  at  Mach 
0.  90  and  at  an  altitude  of  10,  000  feet.  The  resulting  wing  jig  twist  distributions  are 
shown  in  Figure  A-3.  The  final  wing  twist  distribution  for  the  Mach  0.  90  level  flight 
condition  is  presented  in  Figure  A-4. 
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F-lll  AIRPLANE:  WING  SWEEP  =  50° 


Figure  A-3.  Wing  Jig  Twist  Distribution  for  the  F-lll  Composite  Wings 


F-lll  AIRPLANE:  M  =  0.  90,  q  =  824  PSF,  ALT  =  10,000  FT 
WING  SWEEP  =  50° 


Figure  A-4.  Wing  Twist  Distribution  for  the  F-lll  Airplane  with  a 
Composite  Wing  in  Level  Flight  at  Mach  0.  90 
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Figure  A- 5  illustrates  the  shear  load  distribution  on  the  composite  wings  for  the  elastic 
axis  locations  of  20,  44,  and  64  percent  of  the  wing  chord  at  load  factors  of  1.  0  and  5.0. 
Movement  of  the  elastic  axis  as  the  crossply  orientation  is  changed  will  cause  the  span- 
wise  distribution  of  shear  to  increase  with  increasing  load  factor.  The  shear  load  dis¬ 
tribution  for  the  three  elastic  axis  locations  at  the  level  flight  Mach  0.  90  condition  does 
not  change,  because  the  sum  of  the  elastic  twist  and  the  wing  jig  twist  is  identical  for 
each  elastic  axis  location  as  shown  in  Figure  A-4. 

The  effect  of  the  elastic  axis  location  on  the  wing  bending  moment  and  pitching  torque 
perpendicular  and  parallel  to  the  elastic  axis,  respectively,  is  presented  in  Figure 
A- 6. 


F- 111  AIRPLANE:  M  =  0.90s  q  =  824  PSF,  ALTITUDE  =  10, 000  FT 
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Figure  A-5.  Shear  Load  Distribution  for  the  F-lll  Airplane  with  a 
Composite  Wing  at  Load  Factors  of  1.  0  and  5.  0 
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F-lll  AIRPLANE:  M  =  0.90,  q  =  824  PSF,  ALT  =  10,000  FT 


WING  SWEEP  =50° 


Figure  A-6.  Composite  Wing  Bending  Moment  and  Pitching  Torque 
for  the  F-lll  Wing  at  Load  Factors  of  1.  0  and  5.  0 
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The  distribution  of  span  loading  over  the  wing  is  presented  in  Figure  A-7  for  Mach 
0.90,  10,000  feet  and  load  factors  of  1.0  and  5.0.  For  the  level  flight  condition,  the 
span  load  is  identical  for  the  wings  with  the  elastic  axis  located  at  the  20,  44,  and  64 
percent  chordline.  For  increased  load  factors,  the  aft  movement  of  the  elastic  axis 
leads  to  slightly  higher  loading  conditions  near  the  tip. 

The  wing  bending  moment  of  inertia  and  the  wing  polar  moment  of  inertia  are  presented 
in  Figure  A-8  as  a  function  of  the  wing  semispan.  The  ratio  of  bending  moment  of  in¬ 
ertia  to  the  polar  moment  of  inertia  was  0.  60  along  the  wing  semispan. 


F- 111  AIRPLANE:  M  =  0.  90,  q  =  824  PSF,  ALT=10,000  FT 
WING  SWEEP  =  50° 


SEMISPAN 


Figure  A-7.  Span  Load  Distribution  Over  the  F-lll  Wing  with  a 
Composite  Wing  at  Load  Factors  of  1.  0  and  5.  0 
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Figure  A-8.  Wing  Bending  Moment  of  Inertia  and  Polar  Moment 
of  Inertia  for  the  F-lll  Composite  Wings 
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The  elastic  wing  twist  for  the  composite  wing  becomes  more  positive  as  the  elastic 
axis  is  moved  aft.  Wing  section  pitching  moment  does  not  change  with  the  location  of 
the  elastic  axis;  however,  the  aerodynamic  loading  due  to  angle  of  attack  induces  a 
positive  elastic  twist  increment  as  the  elastic  axis  moves  aft,  as  shown  in  Figure  A-9. 

Figure  A-10  illustrates  the  comparison  of  an  aluminum  and  graphite /epoxy  wing  elastic 
twist  distribution  for  the  F-lll  airplane  at  Mach  0.  90.  The  elastic  axis  was  located  at 
the  44%  chordline  of  the  wings.  For  a  constant  load  factor,  the  effect  of  increasing  the 
El/GJ  ratio  of  the  wing  structure  leads  to  reduced  wing  twist  due  to  the  increased 
bending  strength  of  the  wing. 

Vertical  deflection  of  the  composite  wing  is  affected  by  the  location  of  the  elastic  axis. 
Figure  A-ll  illustrates  the  nondimensional  vertical  deflection  of  the  three  composite 
wings  at  load  factors  of  1.  0  and  5.  0  for  a  flight  condition  of  Mach  0. 90.  The  wing  de¬ 
flection  measured  along  the  elastic  axis,  which  has  been  nondimensionalized  with 
respect  to  the  wing  mean  aerodynamic  chord,  decreases  with  an  aft  movement  of  the 
elastic  axis. 

A.  2  LONGITUDINAL  AEROELASTIC  AERODYNAMIC  CHARACTERISTICS  OF  THE 
F-lll  AIRPLANE  WITH  A  COMPOSITE  WING 

The  aerodynamic  methods  used  to  evaluate  the  aerodynamic  consequences  of  the  aero- 
elastic  differences  between  the  airplane  with  a  composite  wing  and  the  solid  steel  IHRT 
model  are  the  same  as  used  in  the  preceding  analysis  of  this  report.  The  assumption 
of  a  linear  twist  distribution  was  made  to  be  consistent  with  the  aerodynamic  methods. 

For  this  effort,  a  composite  F-lll  wing  has  been  analyzed  and  compared  with  the  solid 
steel  HIRT  model  wing.  Three  locations  of  the  elastic  axis  were  evaluated.  Figure 
A-12  illustrates  the  effect  of  load  factor  on  wing  tip  twist  for  the  composite  wings,  the 
aluminum  wing,  and  the  steel  model  wing.  For  the  composite  wing  with  the  elastic  axis 
located  at  the  20%  chordline,  the  tip  twist  is  almost  identical  to  the  aluminum  wing. 

For  the  44%  chord  and  64%  chord  locations  of  the  elastic  axis,  the  maneuvering  elastic 
twist  is  much  smaller  in  magnitude  than  that  of  the  aluminum  wing.  For  the  64%  chord 
elastic  axis  location,  the  maneuvering  elastic  twist  is  positive  with  increasing  load 
factor;  this  can  also  be  seen  from  Figure  A-10.  From  Figure  A-12,  the  wing  tip  twist 
for  a  load  factor  of  5.  0  results  in  a  difference  of  about  3.  5  degrees  between  the  com¬ 
posite  wing  under  flight  conditions  at  Mach  0.  90  and  the  solid  steel  wing  of  the  HIRT 
model  under  equivalent  Reynolds  number  testing. 

For  a  composite  wing  with  the  elastic  axis  located  at  the  20%  chordline,  the  differential 
tip  twist  between  the  model  and  the  airplane  results  in  the  aerodynamic  differences 
presented  in  Figure  A-13,  At  a  one-g  flight  condition,  the  model  wing  and  airplane 
were  matched  at  a  lift  coefficient  of  about  0. 16.  Since  the  model  wing,  under  equiva¬ 
lent  Reynolds  number  testing,  twists  more  under  the  increased  load  factors  than  the 
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Figure  A-9.  Elastic  Wing  Twist  for  the  F-lll  Composite 
Wings  at  Load  Factors  of  1.0  and  5.0 
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WING  SWEEP  *  50° 

F-lll  AIRPLANE  WITH  COMPOSITE  WING  AT  M=0.  90, 
q  =  824  PSF,  ALTITUDE  =  10,  000  FT. 

F-lll  HIRT  MODEL,  M-0.90,  q  =  7891  PSF,  TT  =  300° K 
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Figure  A- 10.  Comparison  of  Elastic  Twist  Distribution  for  a  Composite  and  Aluminum 
Wing  on  the  F-lll  Airplane  at  Load  Factors  of  1.  0  and  5.  0 
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WING  SWEEP  =  50° 

F-lll  AIRPLANE  WITH  COMPOSITE  WING  AT  M*0.90, 
q  =  824  PSF,  ALTITUDE  b  10,000  FT. 

F-lll  HIRT  MODEL,  M  -  0.90,  q-  7891  PSF,  TT  =  300°  K 


Figure  A-ll.  Vertical  Deflection  for  the  F-lll  Composite 
Wings  at  Load  Factors  of  1.  0  and  5.  0 
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WING  SWEEP  =  50° 

F-lll  AIRPLANE  WITH  COMPOSITE  WING  AT  M  =  0.  90,  q  =  824  PSF 
ALTITUDE  =  10,000  FT 
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Figure  A-12.  Effect  of  Load  Factor  on  the  F-lll  Tip  Twist  at 
Mach  0.  90  and  Full  Scale  Reynolds  Number 

composite  wing  at  the  equivalent  flight  condition,  the  model  test  data  would  indicate  a 
lower  lift,  drag,  and  pitching  moment  by  the  amount  shown  in  Figure  A-13.  These 
differences  are  considered  inconsequential  and  are  probably  within  the  balance  accuracy 
band.  Figure  A- 14  shows  the  same  aerodynamic  increments  as  a  function  of  the  angle 
of  attack  for  the  composite  wing  designed  with  the  elastic  axis  at  20%  chord. 


Incremental  aerodynamic  coefficients  are  presented  in  Figure  A-15  for  the  composite 
wing  with  the  elastic  axis  located  at  the  44%  chord  location.  From  Figure  A-12,  the 
differential  tip  twist  was  about  2.2  degrees.  At  a  load  factor  of  5.  0,  the  lift  increment 
due  to  the  differential  tip  twist  is  about  0.0120,  with  the  drag  increment  increased  to 
0.0032.  The  aerodynamic  increments  as  a  function  of  angle  of  attack  are  presented  in 
Figure  A-16. 

Figure  A-17  illustrates  the  incremental  aerodynamic  coefficients  for  the  condition  of 
the  composite  wing  with  the  elastic  axis  located  at  the  64%  chord.  The  differential  tip 
twist  of  3.  5  degrees  at  a  load  factor  of  5.  0  results  in  a  lift  increment  of  0.  0220  and  a 
drag  increment  of  0.0054.  Figure  A-18  illustrates  the  aerodynamic  characteristics 
as  a  function  of  angle  of  attack. 
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WING  SWEEP  =  50° 

F-lll  AIRPLANE  WITH  COMPOSITE  WING  AT  M  =  0.  90,  q  =  824  PSF, 
ALTITUDE  =  10,000  FT 

F-lll  HffiT  MODEL,  M  =  0,  90,  q  =  7891  PSF,  Tt  =  300°K 
SOLID  STEEL  WING.  FULL.  SCALE  REYNOLDS  NUMBER 


LOAD  FACTOR  (g's) 


Figure  A-13.  Estimated  Effect  of  Load  Factor  on  the  Incremental  Aerodynamic  Data 
for  the  F-lll  Study  Configuration  at  Mach  0.90  and  Full  Scale  Reynolds 
Number  (Altitude  =  10,  000  ft)  (Composite  Airplane  Wing  with  Elastic 
Axis  at  20%  Chord) 
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WING  SWEEP  =  50° 

F-lll  AIRPLANE  M  =  0.  90,  q  =  824  PSF,  ALT=10,000  FT 
COMPOSITE  WING  WITH  ELASTIC  AXIS  AT  20%  CHORD 
F-lll  HIRT  MODEL  M  =  0.  90,  q=  7891  PSF,  TT  =  300°K 
SOLID  STEEL  WING.  FULL  SCALE  REYNOLDS  NUMBER 
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Figure  A- 14.  Estimated  Effect  of  Angle  of  Attack  on  the  Incremental  Aerodynamic 
Data  for  the  F-lll  Study  Configuration  at  Mach  0.  90  and  Full  Scale 
Reynolds  Number  (Altitude  =  10,000  ft)  (Composite  Airplane  Wing  With 
Elastic  Axis  at  20%  Chord) 
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WING  SWEEP  =  50° 

F-lll  AIRPLANE  M  =  0.90,  q  =  824  PSF,  ALT  =  10,000  FT 
COMPOSITE  WING  WITH  ELASTIC  AXIS  AT  44%  CHORD 
F-lll  HIRT  MODEL  M  =  0.90,  q=7891  PSF,  TT  =  300°K 
SOLID  STEEL  WING,  FULL  SCALE  REYNOLDS  NUMBER 
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Figure  A-15.  Estimated  Effect  of  Load  Factor  on  the  Incremental  Aerodynamic  Data 
for  the  F-lll  Study  Configuration  at  Mach  0.  90  and  Full  Scale  Reynolds 
Number  (Altitude  =  10,  000  ft)  (Composite  Airplane  Wing  with  Elastic 
Axis  at  44%  Chord) 
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WING  SWEEP  =  50° 

F-lll  AIRPLANE  M  =  0.  90,  q  =  824  PSF,  ALT  =  10, 000  FT 
COMPOSITE  WING  WITH  ELASTIC  AXIS  AT  44%  CHORD 
F-lll  HERT  MODEL  M  =  0.  90,  q  =7891  PSF,  TT  =  300°K 

SOLID  STEEL  WING.  FULL  SCALE  REYNOLDS  NUMBER 
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Figure  A-16.  Estimated  Effect  of  Angle  of  Attack  on  the  Incremental  Aerodynamic 
Data  for  the  F-lll  Study  Configuration  at  Mach  0.  90  and  Full  Scale 
Reynolds  Number  (Altitude  =  10,000  ft)  (Composite  Airplane  Wing  With 
Elastic  Axis  at  44%  Chord) 
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WING  SWEEP  =  50° 

F-lll  AIRPLANE  M  =  0.90,  q  =  824  PSF,  ALTITUDE  =  10,000  FT 
COMPOSITE  WING  WITH  ELASTIC  AXIS  AT  64%  CHORD 
F-lll  HIRT  MODEL,  M  =  0.  90,  q  =  7891  PSF,  TT  =  300°K 
SOLID  STEEL  WING.  FULL  SCALE  REYNOLDS  NUMBER 


-0.0500 
-0.  0400 
ACm  -0.0300 
-0. 0200 
-0.0100 
0 


1.0  2.0  3.0  4.0  5.0  6.0 

LOAD  FACTOR  (g's) 


Figure  A-17.  Estimated  Effect  of  Load  Factor  on  the  Incremental  Aerodynamic  Data 
for  the  F-lll  Study  Configuration  at  Mach  0.  90  and  Full  Scale  Reynolds 
Number  (Altitude  =  10,000  ft)  (Composite  Airplane  Wing  with  Elastic 
Axis  at  64%  Chord) 


/ 


121 


AEDC-TR-75-62 


WING  SWEEP  =  50° 

F-lll  AIRPLANE  M  =  0.  90,  q  =  824  PSF,  ALT  =  10,000  FT 
COMPOSITE  WING  WITH  ELASTIC  AXIS  AT  64%  CHORD 
F-lll  HIRT  MODEL,  M  =  0.  90,  q  =  7891  PSF,  TT  =  300°K 
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Figure  A-18.  Estimated  Effect  of  Angle  of  Attack  on  the  Incremental  Aerodynamic 
Data  for  the  F-lll  Study  Configuration  at  Mach  0.90  and  Full  Scale 
Reynolds  Number  (Altitude  =  10,000  ft)  (Composite  Airplane  Wing  With 
Elastic  Axis  at  64%  Chord) 
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A3.  ANALYTICAL  RESULTS 

The  results  of  this  analysis  indicate  that  the  simulation  of  a  composite  wing  by  means 
of  a  solid  steel  model  wing  for  equivalent  Reynolds  number  testing  will  require  an  addi¬ 
tional  wing  or  wings  to  be  jig  twisted  to  obtain  the  desired  level  of  model/airplane 
similarity  at  a  particular  flight  condition.  This  will  of  course  depend  on  the  importance 
of  the  maneuvering  aerodynamic  characteristics  for  various  concepts  tested,  as  well 
as  the  importance  of  maintaining  the  similarity  in  Reynolds  number. 

From  this  analysis,  the  wing  stiffness  ratio  (El/GJ)  and  the  location  of  the  elastic 
axis  simulated  by  the  solid  steel  model  wing  will  influence  the  number  of  wings  re¬ 
quired  to  simulate  the  airplane  through  the  performance  envelope.  For  a  composite 
wing  with  the  elastic  axis  located  near  the  forward  portion  of  the  airfoil  section,  a 
single  wing  may  be  jig  twisted  to  provide  good  test  data  through  an  angle -of-attack 
range  that  would  be  equivalent  to  load  factors  up  to  5.  0  for  a  combat  airplane.  As  the 
elastic  axis  is  moved  aft,  the  simulation  by  a  solid  steel  model  wing  becomes  more 
unacceptable.  The  incremental  aeroelastic  deformations  between  the  solid  steel  wing 
and  the  composite  wing  result  in  moderate  effects  on  the  aerodynamic  data. 

The  elastic  axis  of  future  composite  structures  should  move  aft  of  the  isotropic 
position  to  achieve  the  best  compromise  of  flutter  and  divergence  speeds.  The  flutter 
speed  will  increase  with  the  aft  movement  of  the  elastic  axis,  all  other  flutter  variables 
remaining'  unchanged,  and  the  divergence  speed  will  decrease.  Current  isotropic 
designs  are  rarely  critical  in  divergence  when  a  satisfactory  flutter  margin  over  the 
flight  envelope  has  been  attained. 

The  optimum  amount  of  aft  movement  of  the  elastic  axis  can  only  be  determined  by 
examining  both  the  flutter  and  divergence  speeds  for  two  or  more  elastic  axis  loca¬ 
tions.  The  aft  position  (64%C)  considered  in  this  work  was  selected  for  the  purpose 
of  studying  the  effects  on  the  aerodynamic  response  and  represents  an  extreme  position. 
It  may  well  be  that  a  practical  limit  of  shifting  the  elastic  axis  aft  will  be  between  44 
and  64  percent  of  the  local  chord.  It  is  clear  from  the  results  of  this  report  that,  if 
a  large  aft  movement  of  the  elastic  axis  of  the  full  size  anisotropic  structure  is  needed 
to  reduce  wing  weight,  one  or  more  additional  model  wings  could  be  required  to  obtain 
the  desired  level  of  model/airplane  simularity. 


From  the  previous  study,  several  testing  techniques  may  be  used  to  improve  the  simu¬ 
lation  of  the  composite  wing  with  a  steel  model  wing.  Tunnel  temperature  reduction 
will  allow  the  full-scale  Reynolds  number  to  be  simulated  with  a  reduced  dynamic  pres¬ 
sure.  This  would  result  in  reducing  the  model  wing  deflection  at  a  selected  equivalent 
flight  condition.  Mismatching  Reynolds  number  would  be  the  second  method,  which 
again  results  in  reducing  the  tunnel  dynamic  pressure  and  thus  the  model  wing 
deflections. 
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The  composite  wing  deflections  that  have  been  analyzed  during,  this  study  were  based 
on  a  wing  stiffness  ratio  of  3.0.  For  conditions  with  reduced  wing  stiffness  ratios,  the 
wing  deformations  will  begin  to  approach  the  aluminum  wing  characteristics.  The  in¬ 
creased  wing  bending  will  lead  to  an  increased  wing  twist.  This  condition  will  improve 
the  ability  to  maintain  model  wing  similarity  with  the  flight  vehicle  at  the  equivalent 
Reynolds  number  condition. 

Chordwise  bending  of  the  wing  structure  has  not  been  considered  in  the  study.  The 
aerodynamic  analysis  is  based  on  the  twist  distribution  of  the  wing  only.  As  stated 
previously,  the  wing  bending  characteristics  are  an  integral  part  of  the  wing  twist  dis¬ 
tribution,  thus  the  aerodynamic  characteristics  resulting  from  a  nonplanar  wing  are  not 
considered  directly. 

A4,  CONCLUSIONS 

This  study  indicates  that,  for  equivalent  Reynolds  number  testing  in  HIRT,  multiple 
steel  wings  may  be  required  to  simulate  the  composite  wing  shape  accurately  and  thus 
obtain  good  experimental  data.  Several  considerations  qualify  this  result.  First,  when 
there  is  a  close  agreement  between  the  model  wing  and  the  airplane  composite  wing 
elastic  axis  location,  appreciable  twist  differences  occur,  producing  aerodynamic  data 
that  does  not  accurately  represent  the  characteristics  of  the  aircraft.  For  the  subject 
configuration,  the  composite  wing  had  a  stiffness  ratio  (El/GJ)  of  3.0,  which  is  high 
relative  to  an  aluminum  structure  of  about  1.1.  As  the  composite  wing  stiffness  ratio 
is  reduced  to  levels  of  the  aluminum  structure,  one  model  wing  may  be  sufficient  to 
maintain  wing  shape  similarity  throughout  the  pitch  range.  For  the  condition  where 
appreciable  twist  differences  occur,  drag  is  most  affected,  followed  by  pitching  moment 
and  lift.  A  second  consideration  is  the  importance  of  high  g  maneuver  performance. 

For  a  transport  design  at  load  factors  of  2.5,  one  model  wing  may  be  sufficient.  Military 
maneuvering  designs  have  emphasis  on  sustained  maneuver  performance  (L/D)  and  in¬ 
stantaneous  maneuverability  (CLmax,  buffet  onset,  controlability).  This  type  of  design 
is  the  most  likely  to  require  multiple  wing  testing. 

The  model  deformation  can  be  traded  against  Reynolds  number  to  get  model  wing  shape 
similarity  for  nonsimilar  test  Reynolds  numbers  by  constraining  the  test  dynamic  pres¬ 
sure  as  a  function  of  model  wing  shape  rather  than  Reynolds  number.  Control  capability 
of  tunnel  operating  temperature  would  provide  an  additional  dimension  for  the  control  of 
dynamic  pressure,  model  shape,  and  test  Reynolds  number. 
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NOMENCLATURE 


Symbol 

Description 

Units 

AEAD 

Slope  of  the  wing  elastic  axis 

deg. 

A  LED 

Slope  of  the  wing  leading  edge 

deg. 

ATED 

Slope  of  the  wing  trailing  edge 

deg. 

AR 

Wing  aspect  iatio 

dim. 

A 

W 

Wing  downwash  matrix 

-1 
in  . 

A 

F 

Wing  image-vortex  matrix 

in  . 

BL 

Buttock  line 

in. 

b 

Wing  span 

in. 

C 

Wing  mean  aerodynamics  chord 

in. 

CD 

Drag  coefficient 

dim. 

C 

-1 

DEq? 

Drag  coefficient  of  external  store 

deg  . 

Cdl 

Induced  drag  coefficient 

dim. 

^Cn 

Incremental  drag  coefficient  due  to  airplane/model  wing 

i-/ 

twist  differences 

dim. 

CT 

Lift  coefficient 

dim. 

L 

cTT, 

LEa 

External  store  lift  curve  slope 

deg  . 

C 

LFa 

Fuselage  lift  curve  slope 

deg  . 

°la 

Wing  lift  curve  slope 

deg  . 

aCl 

Incremental  lift  coefficient  due  to  airplane/ model  wing 

twist  differences 

dim. 

C 

i 

Wing  section  lift  coefficient 

dim. 

CM 

Pitching  moment 

dim. 

C 

MFcf 

Fuselage  pitching  moment  curve  slope 

deg  1. 
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NOMENCLATURE  (cont'd) 

Symbol 

Description 

Units 

C 

m 

o 

Wing  section  pitching  moment  at  zero  lift. 

dim. 

AC 

m 

Incremental  pitching  moment  coefficient  due  to  airplane/ 
model  wing  twist  differences 

dim. 

CR 

Wing  root  chord 

in. 

CT 

Wing  tip  chord 

in. 

f 

E 

Modulus  of  elasticity 

Ib/in 

EA 

Elastic  axis  location  (%  chord) 

- 

[ea] 

External  store  matrix 

lb"1. 

ETA  (V) 

Location  of  wing  strip  along  the  Y-axis  (%  semispan) 

r 

G 

Shear  modulus  of  elasticity 

lb/in' 

H 

Wing  strip  halfwidth 

in. 

I 

Wing  section  moment  of  inertia 

.  4 
in  . 

H 

Moment  of  inertia  matrix 

.  4 
in  . 

IYCG 

Moment  of  inertia  of  the  airplane  about  the  pitch  axis 
through  the  center  of  gravity. 

in-lb- 

J 

Wing  section  polar  moment  of  inertia 

,  4 

m  . 

a 

Span  load 

lb. 

M 

Mach  number 

dim. 

MS 

Model  station 

in. 

m 

Bending  moment  per  unit  pitching  moment 

dim. 

m 

o 

Lift  curve  slope  including  compressibility  effects 

deg  1 

n 

Number  of  strips  on  the  wing  semispan 

dim. 

nz 

Load  factor  along  Z  axis 

dim. 

Pm 

T 

Horizontal  tail  load 

lb. 

Q(q) 

Free  stream  dynamic  pressure 

lb/ft2 

2 

sec  . 
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NOMENCLATURE  (cont'd) 


Symbol 

Description 

Units 

R 

e 

Reynolds  number 

dim. 

S 

Wing  area 

2 

in  . 

H 

Wing  elasticity  matrix 

.  -1 
in  . 

T 

T 

Tunnel  storage  temperature 

°  K 

t 

Torsional  moment  around  the  elastic  axis  per  unit  pitching 
moment 

dim. 

V 

Free  stream  velocity 

ft/sec. 

K/Vz] 

Overspeed  effect  matrix 

dim. 

L  J 

w 

Airplane  gross  weight 

lb. 

WL 

Waterline 

in. 

w 

p 

Induced  velocity  normal  to  the  wing  chord  plane 

ft /sec. 

XC4 

Fuselage  location  of  the  25%  station  of  the  wing  mean 
aerodynamic  chord 

in. 

YC4 

Buttock  line  location  of  the  25%  station  of  the  wing  mean 
aerodynamic  chord 

in. 

ZC4 

Waterline  location  of  the  25%  station  of  the  wing  mean 
aerodynamic  chord 

in. 

ds 

Differential  distance  along  the  elastic  axis 

in. 

de_ 

da 

Downwash  at  the  horizontal  stabilizer 

dim. 

a 

Angle  of  attack  of  the  airplane 

deg. 

aF 

Angle  of  attack  of  a  wing  section  including  aeroelastic  effects 

deg. 

a. 

inc 

Incidence  angle 

deg. 

a  T 
oL 

Wing  section  zero  lift  line  angle  of  attack 

deg. 

O' 

r 

Wing  root  section  angle  of  attack 

deg. 

a 

Elastic  angle  of  attack  due  to  bending  and  torsional 

s 

moments  along  the  elastic  axis 

deg. 
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NOMENCLATURE  (cont'd) 


Symbol 

Description 

Units 

Of 

r 

Wing  jig  twist  relative  to  the  root  chord 

deg. 

r 

Circulation 

ft'Vsec, 

6fl 

Flap  deflection 

deg. 

6sl 

Slat  deflection 

deg. 

6sp 

Spoiler  deflection 

deg. 

®CG 

Pitch  acceleration  of  airplane 

deg/sec2. 

A 

Wing  sweep  angle 

deg. 

P 

Fluid  density 

lb  sec2/ft 
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